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1 Introduction 



In my Les-Houches lectures of 1982 I described the inverse scattering method of 
solving the integrable field-theoretical models in 1 + 1 dimensional space-time. Both 
classical case, stemming from the famous paper by Gardner, Green, Kruskal and 
Miura of 1967 on KdV equation, and its quantum counterpart, developed mostly by 
Leningrad group around 1978-79, were discussed. In particular, the algebraic way of 
deriving the Bethe-Ansatz equations was presented, but its use was not illustrated 
in any detail. I just stated in the end of the lectures, that "The hard work just 
begins" . 

In this course I shall exactly describe this work. During last 10 years I lectured 
several times on this subject, but this particular lecture course is the longest and 
more detailed. 

In the announcement of the School my course is called "Hidden symmetries in 
integrable models". The term "symmetry" in modern literature on mathematical 
physics is supplied with several adjectives, such as hidden, dynamical, broken, de- 
formed etc., but there is no exact definition for all this. Thus I decided to change my 
title in the proceedings to reflect more adequately the actual content of the course. 

The term "integrable models" in the title refers to particular family of quantum 
field-theoretical models in 1 + 1 dimensional space-time, which are solvable by means 
of the quantum variant of inverse scattering method. The adjective "integrable" 
stems from a paper of Zakharov and me of 1971 where the KdV equation was 
shown to allow the interpretation as an integrable (though infinite-dimensional) 
hamiltonian system. 

The most famous example of the integrable model is the Sine-Gordon equation 

□y? + m 2 sin (p = (1) 

for the scalar field ip(x,t), which is relativistic and nonlinear. 

As is well known, the quantum field-theoretic models with interaction are plagued 
by infinities and some regularization is necessary. The discretizing of the space, or 
going to the lattice, is one way for such a regularization. It reduces the field model 
in the finite volume to a system with finite number of degrees of freedom. 

In the beginning of 80-ties, due mainly to the work of Izergin and Korepin it was 
realized, that the integrable models allow the lattice counterparts, which are also 
integrable and can be interpreted as the quantum spin models of magnetic chains. 
This unexpected but very welcome connection showed the universality of the spin 
chains in the domain of integrable models. In my course I shall begin with and 
speak in detail on the spin chains. The field-theoretical models will appear as their 
particular continuous space limits. 

One can ask, what is good in 1 + 1 models, when our space-time is 3 + 1- 
dimensional. There are several particular answers to this question. 

1. The toy models in 1 + 1 dimension can teach us about the realistic field- 
theoretical models in a nonperturbative way. Indeed such phenomena as renor- 
malization, asymptotic freedom, dimensional transmutation (i.e. the appearance 
of mass via the regularization parameters) hold in integrable models and can be 
described exactly. 



2. There are numerous physical applications of the 1 + 1 dimensional models in 
the condensed matter physics. 

3. The formalism of integrable models showed several times to be useful in 
the modern string theory, in which the world sheet is 2-dimensional anyhow. In 
particular the conformal field theory models are special massless limits of integrable 
models. 

4. The theory of integrable models teaches us about new phenomena, which were 
not appreciated in the previous developments of Quantum Field Theory, especially 
in connection with the mass spectrum. 

5. I cannot help mentioning that working with the integrable models is a de- 
lightful pastime. They proved also to be very successful tool for the educational 
purposes. 

These reasons were sufficient for me to continue to work in this domain for last 
25 years (including 10 years of classical solitonic systems) and teach quite a few of 
followers, often referred to as Leningrad-St. Petersburg school. 

I am very grateful to the organizers of the school Professors A.Connes and 
K.Gawedsky for inviting me to lecture. First, it is always nice to be in Les-Houches 
and it is already my third lecture course here. Second, the texts of two previous lec- 
tures were transformed into monographs later. I hope that this course will eventually 
lead to one more book dedicated to the quantum theory of solitons. 

2 General outline of the course. 

My usual methodological trick in teaching is not to begin in full generality but rather 
to choose a representative example and explain on it all technical features in such a 
way, that generalization become reasonably evident. Thus I begin the course with 
the concrete example of the magnetic model - the spin 1/2 XXX chain. 

All the ingredients of Algebraic Bethe Ansatz, which is another name for the 
Quantum Inverse Scattering method, namely Lax operator, derivation of Bethe- 
Ansatz equations, thermodynamic limit, will be presented in full detail on this ex- 
ample. Then the generalization to XXX model of spin s > 1/2 and XXZ model 
will follow. Here only features, which distinguish this model from the basic one, 
will be described. Finally the continuous field-theoretical models such as Nonlinear 
Schroedinger Equation, S 2 nonlinear a-model, principal chiral field model for SU (2) 
and Sine-Gordon model will be included as particular limits of spin chains. 

We will see, that in the description of dynamics the finite time shift 

U = e~ tHA (2) 

will appear naturally. This will make our discretization scheme more consistent, 
time and space being discrete simultaneously. 

Now I shall present some kinematics of the models we shall consider. 

As "space" we shall consider a discrete circle, namely the ordered set of points, 
labeled by integers n with the identification n = n + N, where N is a fixed positive 
integer. As a "fundamental domain" we shall take n = 1,...,N. The integer N 
plays the role of the volume of the space; the identification reflects the periodic 
boundary condition. 



Formal continuous limit will be described by introducing the lattice spacing A 
and coordinate x = nA which becomes continuous in the limit A — > 0, N — > oo. In 
particular the following rule will be used 

% = *(*-!/) , (3) 

so that the Kroneker symbol 5 mn is of order A. 

The quantum algebra of observables A is generated by dynamical variables I™, 
attached to each lattice site n. Index a assumes some finite number of values. The 
algebra A is defined by fixing the set of commutation relations between X". These 
relations are called ultralocal, when X^ and X@ commute for n ^ m. A more 
relaxed condition of locality is that X^ and X@ do not commute only for \n — m\ 
small, in particular for n = m — 1 and n = m + 1 . 

Let us present examples, beginning with ultralocal case. 

1. Canonical variables ip°, 7r°, a — 1, . . . , I with the relations 

[<ftM = o, K>^] = o, (4) 

Ki ^n] = ihISnmSaf3 ■ (5) 

Here [ , ] is used for commutator 

[a, b] = ab — ba , (6) 

h is a Planck constant, which we soon shall drop, I is as unity in algebra A. 
We can call / a number of degrees of freedom per lattice site; the full number 
of degrees of freedom is Nl. 

Each pair of canonical variables is represented in an infinite dimensional Hilbert 
space, which can be chosen as Z/2(R) with tp and ir being operators of multi- 
plication and differentiation 

v /(0 = C/(0 ; vr /(0 = -^/(0 • (7) 



2. Spin variables S 1 ", a — 1, 2, 3 with the relations 

[S m , S n ] = ilhe a fs 7 S^S mn , (8) 

where e al 3 y is a completely antisymmetric tensor, £123 = 1. 

Mathematically these variables define a Lie algebra s/(2), the finite dimensional 
representations of which are labeled by half-integer s — 0, 1/2, 1, ... and are 
realized in C 2s+1 . In the smallest nontrivial dimension (s = 1/2) operators S% 
are represented by Pauli matrices a a 




as = h/2a a . 



3. Weyl variables. 

For one degree of freedom per lattice site these variables consist of a pair u n , 
v n with the exchange relations 

^m^n ^n^m i ^mVn ^nVm j ^"mVn ^nV"m j 777. 7^ 71 , 

w n t>„ = qv n u n , (10) 

where g is a given complex number. Usually it assumes the values on a circle 
and is parametrized by a real number 7 

q = e lh ~ / . (11) 

4. g-deformed spin variables S%. In writing the commutation relations we shall 
drop the lattice index n and present them for fixed n as follows 

q s3 S ± = q ±l S ± q s3 ; (12) 

[S\S-]= {q ] l\ ] , (13) 

via generators denoted by q s , S + and S~ . When q — > 1 (or 7^0) these 
relations turn into the usual s/(2) relations for S a , S' ± being the usual combi- 
nations 5± = ± iS 2 . 

Thus they define a g-deformed s/(2) algebra, denoted by sl q (2). For generic 
g the finite dimensional representations are given in the same spaces C 2s+1 
as in nondeformed case. However for q on the circle some new interesting 
representations occur, more on this below. 

The Hilbert space for the representation of ultralocal algebra A has a natural 
tensor-product form 

N 

H = Yl ®K = h x <g) h 2 <S> ■ ■ ■ <S> K <S> ■ ■ ■ <S> h N (14) 

n=l 

(where all h n could be the same) and variables X% act nontrivially only in the space 
K 

^ = I®I(8)... (g)X a (g) ...(g)! . (15) 
n-th place 

In the continuous limit we encounter the problem of considering the infinite tensor 
products, which is known to be intricated from the time of v. Neumann. In fact, 
the concrete examples of the continuous limits will give the instances of the rather 
nontrivial constructions of such products. 

The simplest example of the nonultralocal relation is furnished by the exchange 
algebra for variables w n with the relations 



(16) 



The Hilbert space for this algebra is also a tensor product, but with the length 
being around N/2. We shall not discuss it here, however, the example being given 
just for illustration. 

General considerations a-ld Darboux theorem in classical mechanics state that 
the canonical variables are generic in the sense that all other types of dynamical 
variables can be expressed through them. Let us illustrate it on the examples above. 

Let ip, if)* be a complex pair of canonical variables with relation 



[V,V*]=I (17) 
(i.e. ip = A%((p + wr), if)* = ^75 (y 9 — wr))i then the variables 

S+ = S 1 + iS 2 = if)*{2s-if)*iP) ; (18) 

S- = S 1 -iS 2 = if); (19) 

S 3 = if;*if)-s (20) 



satisfy the spin commutation relation. Here s is any complex parameter, which can 
be called spin, as the Casimir C = (S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 has value C = s(s + 1). 
Weyl pair u, v can be realized via canonical variables 7r and tp as follows 

u = e iwQ , v = e^ p , 

where P and Q are any complex numbers and q = e % 

Finally, the deformed spin variables q s3 , S ± can be realized as 

S ± = e ±i*/2(i + m 2 e ±2i^ e ±™/2 . (22) 

q s3 = e iw , (23) 

where m 2 and 7 are parameters. 

Of course the Hilbert space for the canonical variables is infinite dimensional and 
the representations for the spin variables appear as reductions. For example, if 2s 
is integer then the subspace consisting of 2s + 1 states u, ip*ui, . . . (-?/>*) 2s u;, where u> 
is a vector annihilated by ip, ipu = 0, is invariant with respect to action of S , S 3 . 

The Weyl pair u, v in form (|2T| ) has finite dimensional representation if q is a 
root of unity, or when PQ/2n is a rational number. Corresponding reduction exists 
also for the deformed spin variables if 7/2-7? is rational and it gives the so-called 
finite dimensional cyclic representation of sl q {2). 

These reductions evidently reflect some global aspects of the would-to-be quan- 
tum Darboux theorem. 

This completes my short introduction to the kinematics of the models I plan to 
consider. I come to the representative example — spin 1/2 XXX chain. 



(21) 



3 XXX x /2 model. Description. 

The origin of the abbreviation XXX will become clear soon. As I have already said 
this is quantum model, defined in a Hilbert space 

N 

H N = J] ®h n , (24) 

71=1 



where each local space h n is two-dimensional 

h n = C 2 . (25) 

The spin variables S° are acting on each h n as Pauli matrices divided by 2. 
There are several important observables, such as the total spin 

n 

or the total hamiltonian 

TT \ ^ ( na n« 

n — I D n D n+l ~~ 7 

where the periodicity 

na n« 

°n+V — D n 

is to be taken into account. We have 

[H, S a ] = , 

which reflects the sl(2) symmetry of the model. 

The abbreviation XXX is used to stress this invariance which is reflected in the 
fact that all coefficients in front of combination S^S^ +1 in hamiltonian are equal. A 
more general hamiltonian 

h = y: J a s:s: +1 (so) 

a,n 

with parameters J a corresponds to XYZ spin 1/2 model. 

Our problem is to investigate the spectrum of H. Of course for finite N it is just 
a problem about the matrix 2 N x 2^ accessible by computer. However we will be 
interested in the limit N — > oo, when only analytic methods work. 

The core of our approach is a generating object, called Lax operator. This object 
is a rather long shot from the Lax operator of KdV equation but historically this 
name was fixed for any linear operator, entering the auxiliary spectral problem of 
the classical inverse scattering method. 

The definition of the Lax operator involves the local quantum space h n and the 
auxiliary space V, which for the beginning will be also C 2 . Lax operator L n a (A) 
acts in h n <g> V and is given explicitly by the expression 

L n , a (X) = \l n ®l a + i^S%®a a , (31) 

a 

where I n , S% act in h n and I a , a a are unit and Pauli matrices in V = C 2 ; A is a 
complex parameter, usually called the spectral parameter, reminding its role as an 
eigenvalue in the original Lax operator. 

Alternatively L nA (X) can be written as 2 x 2 matrix 

WA)=( A + f x %) , (32) 



(26) 

(27) 
(28) 
(29) 



acting in V with entries being operators in quantum space h n . One more form uses 
the fact that operator 

P = i(I®I + ]Ta a ®(X Q ) (33) 

is a permutation in C 2 (g) C 2 

P a®b = b® a . (34) 

In terms of P n ,a ; which makes sense due to the fact that h n and V are the same C 2 , 
we have 

L n ,a(\) = (\- t -)I n , a + lP n , a . (35) 

Now we establish the main property of Lax operator — the commutation relation 
for its entries. As we have four of them we are to write down 16 relations. Our 
convenient notations allow to write them all in one line. 

Consider two exemplars of Lax operators L niQl (A) and L n)(l2 (/z) with the same 
quantum space and V\ and V 2 serving as corresponding auxiliar spaces. The products 
L n , ai W L na2 (fi) and L nja2 (/i) L nai (A) make sense in a triple tensor product h n <g> 
Vi®V2- We claim, that these two products are similar operators with the intertwiner 
acting only in V\ ® V% and so not containing quantum operators. In other words, 
there exists an operator i? ai a2 (A — \i) in V\ ® V2 such that the following relation is 
true 

R ai ,a 2 (^- A*) L nAl (X) L n , a2 (/f) = L nja2 {n) L n>ai (\) i2o li02 (A - n) . (36) 
The explicit expression for R aiiCL2 (X) is 

(A) = AI 

11,02 ~ 01 ai,tt2 ) 

(37) 

where I ai ,a 2 an d Pai,a 2 are unity and permutation in V\ <8> Vi- 

Comparing Q3"5j ) and (P7|) we see that the Lax operator L n a (A) and operator 
Rai,a 2 (tyi which we shall call i?-matrix, are essentially the same. 

To check (j36|) it is convenient to use the form (|37|) for L n a (A) and the commu- 
tation relation for permutations 

Pn,ai Pn,d2 Pai,a2 Pn.ai Pn,a 2 Pa2,ai (38) 

together with evident symmetry 

Pa2,ai Pai,a 2 ■ (39) 

The importance of the relation ([36]) will become clear momentarily. We shall call 
it in what follows the fundamental commutation relation (FCR). Its general place 
in the family of the Yang-Baxter relations will become evident later. 

The Lax operator L n ^ a (X) has a natural geometric interpretation as a connection 
along our chain, defining the transport between sites n and n+1 via the Lax equation 

ip n+1 = L n ip n (40) 

{ ib 1 \ 

for vector ip n — I ,5 with entries in Ji. The ordered product over all sites 



between n 2 and n\ 



TZaW = Ln 2 ,aW ■ ■ ■ V,a(A) (41) 



defines the transport form n\ to ri2 + 1 and the full product 

T N JX) = L N JX) . . . LUX) (42) 



JV,a 



(43) 



is a monodromy around our circle. The last operator is given as a 2 x 2 matrix in 
the auxiliary space V 

' Ajv(A) , B N (X) \ 
v Cjv(A) , -Dat(A) ) 

with entries being operators in the full quantum space 7i. As in classical case the 
map from local dynamical variables ( S% ) to monodromy T/v, a (A) is a tool for 
solving the dynamical problem. We shall see that T Na (X) is a generating object for 
main observables such as spin and hamiltonian, as well as for the spectrum rising 
operators. 

For that we shall establish first the FCR for T/v, a (A). We claim, that it has 
exactly the same form as the local FCR ([36]), namely 

Ra u a 2 (X - fl)T ai (X)T a . 2 (fi) = T a2 (fi)T ai (X)R aua2 (^ - ■ ( 44 ) 

We dropped here the index N and shall do it in what follows as soon as it does not 
lead to confusion. 

Derivation of ([0]) is very simple and uses the advantages of our notations. We 
shall prove it for any transport operator. It is clear, that it is enough to consider the 
transport along two sites, i.e. n — ■> n + 2. With short notations i? aijCl2 (A — fi) = Ri 2 , 
L n , ai (X) = L\, L n+1>ai (X) = Li, L„, a2 (/i) = L 2 , L n+ iUv) = L 2 we nave 

Ri2L' 1 LiL' 2 L2 = (due to commutativity of L\ and L' 2 ) 

= Ri2L' 1 L' 2 LiL2 = (due to the local FCR for L\, L 2 and L[, L' 2 ) 

= L' 2 L / 1 L 2 Lii?i2 = (due to commutativity of L[ and L 2 ) 
= L 2 L 2 L' 1 LiRi2 

This completes the proof. 

The monodromy Tjv, a (A) is a polynomial in A of order iV 

T N , a (X) = X N + iX N -' Y,(S a ® O + . . . , (45) 

a 

so that total spin S a appears via the coefficient of next to the highest degree. Now 
we shall find the place for the hamiltonian. The FCR ((Pj) shows, that the family of 
operators 

F(X) = trT(A) = A(X) + D(X) (46) 

is commuting 

[F(A),F(a*)] = . (47) 
Its nontrivial A expansion begins with power A^" 2 

N-2 

F(X) =2X N +J2 Qi^ 1 (48) 

1=0 

and produces N — 1 commuting operators Qi. We shall show, that H belongs to 
this family. 



For this note, that the point A = i/2 is rather special 

L n . a {i/2) = zP„, a (49) 

and of course for any A 

-jrL na (X) = l n , a . (50) 
aX 

This makes it easy to control the expansion of F(\) in the vicinity of A = i/2. 
We have 

T N;a (i/2) = i N P Nja Pjv-l,a • • • Pl,a • (51) 

This string of permutations is easily transformed into 

Pl,2 P2,3 • • • ~Pn-1,N PjV,a 

(52) 

by taking permutations one after another from left to right and taking into account 
the properties ( |3"%D and and commutativity of permutations with completely 
different indeces. Now the trace over the auxiliary space is easily taken 

tr Pjv, a — IiV j (53) 

so that 

U = r N tv a T N (i/2) = P lj2 P 2>3 . . . Pat-w • (54) 

It is easy to see, that U is a shift operator in TC. The property (0) can be rewritten 

as 

Pni,ri2 -X-ri2 Pni,ri2 -^-ni ■ (55) 

Thus 

X n U = Pi_2 • • • X n P n _i jn P n ,n+1 • • • PjV-l,AT = 
Pl,2 • • • Pn-l,n X n -l P n ,n+1 • • • PjV-l,JV = UX n ^ . (56) 

Operator U is unitary 

U*U = UU* = I , (57) 
because the permutations have properties 

P* = P ; P 2 = I , (58) 

and we have 

U- x X n U = X n _! , (59) 

which allows to introduce important observable — momentum. By definition, mo- 
mentum P produces an infinitesimal shift and on the lattice it is substituted by shift 
along one site 

e iP = U . (60) 
We proceed now to expand F(\) in the vicinity of A = i/2. We get 



= ^" 1 ^P iVia ...P n , a ...P lia , (61) 

\=i/2 n 



where ~ means that corresponding factor is absent. Repeating our trick we transform 
this after taking the trace over V into 



dX 



Fa(X) 



;N-1 



X=i/2 



^ Pl,2 • • • Pn-l,n+l • • • P 



N-1,N 



(62) 



We can cancel most of permutations here, multiplying by U as a result we get 



d 

dX 



F a {X) F^X)- 1 



X=i/2 



d 

dX 



lnF a (A) 



X=i/2 



71+1 



(63) 



Using ( |33"D we can rewrite the expression (|27| ) for the hamiltonian as 

H 



2 n,n+l 2 



(64) 



and comparing (S3) and ( |63| ) we have 



2dA 



lnF(A) 



X=i/2 



N 
~2 



(65) 



Thus I have shown, that if indeed belongs to the family of AT — 1 commuting 
operators generated by the trace of monodromy F(X). One component of spin, 
say S 3 , completes this family to N commuting operators. This is a proof of the 
integrability of the classical counterpart of our model, which can be considered as a 
system with N degrees of freedom. In the next section I shall show, how to describe 
the set of eigenstates of commuting family F(X) using the offdiagonal elements of 
monodromy. 



4 model. Bethe Ansatz equations. 

Here I shall describe a procedure to diagonalize the whole family of operators F(X) 
in a rather algebraic fashion, based on the global FCR (f44|) and some simple proper- 
ties of local Lax operators. In a way the working generalizes the simplest quantum 
mechanical treatment of harmonic oscillator hamiltonian n = %p*%p based on com- 
mutation relations = I and existence of a state uj such that ipu = 0. 
Let us write the relevant set of FCR 

[5(A), £(//)] = ; (66) 
A(X) Biji) = f(X - /i) B{p) A(X) + g(X - /*) B(X) AQi) ; (67) 
D(X) B(fj) = h(X - y) B(fj) D(X) + k(X - //) B(X) D{ji) , (68) 



where 



/(A) = ^ ; 9(A) = { ; 

h(X) = ^±1 ; k(X) = - l - . (69) 



To read these relations from one line FCR fl44|) one must use the explicit matrix 
representation of FCR in V <g> V. All objects in it are 4x4 matrices in a natural 
basis 

e 1 = e + <S> e + , e 2 = e + £g> e_ , e 3 = e_ <S> e + , e 4 = e_ <8> e_ 
in C 2 ® C 2 , where 



1 





1 



Matrix P assumes the form 



/ 1 
















1 








1 








I o 








1 J 



so that R(X) looks like 



R(X) 



( a(X) 



V 



b(X) c(A) 
c(A) 6(A) 



a(A) / 



(we do not write in the zeros), where 

a = A + % , b = X , c 
The matrices T ai (A) and T a2 (/i) take the form 



r ai (A) 



/ A(A) 5(A) 

A(A) B(X) 

C(X) D(X) 

C(X) D(X) J 



and 



MO*) 5 M 
c(a*) £>(/*) 

V 



\ 

C(//) D(/x) / 



Thus we have 



T ai (X)T a2 (fx) 



( A{\)A(fi) A(X)B(ji) B(X)A^) B{X)B{^) \ 

A(X)C(fi) A(X)D(n) B(X)C(fi) B(X)D(fi) 

C{X)A{fi) C(X)B(n) D(X)A(n) D(X)B(n) 

V C(A)C(//) C{X)D{fj) D(X)C(fx) D(X)D(fj) j 



(70) 
(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



and T a2 (n) T ai (X) is given by matrix, where in all matrix elements the factors have 
opposite order. Now to get ( |6"T|) one is to use the (1,3) relation in FCR ([4*4]) 



a(A - /i) B(X) A(fj) = c(A - //) £(//) A(A) + 6(A - //) 5(A) (78) 



and interchange A <-> fi to get 

Other relations are obtained similarly. 

The exchange relations (|6Tf) and (|B8l ) substitute the relations 

^ n = ( n + l)V>, = (n-l)-0* (80) 

for the harmonic oscillator with n = Now the analogue of "highest weight" u 

such as i/j uj = will be played by a reference state Q such, that 

C{X)Q = . (81) 

To find this state we observe that in each h n there exists a vector u n such that the 
Lax operator L n>a (X) becomes triangular in the auxiliary space, when applied to it 

L n (X)u n = ( X + 5 x *_^ oo n . (82) 

This vector is given by u n = e + . By * we denote operator expressions which are not 
relevant for us. Now for vector Q in H 

= T[®uj n (83) 

n 

we get 



a N (X) 



* 



r(A)n = ^ W ^ (A) j n , (84) 

where 

«(A)=A + ^, 5(A) =A-i. (85) 

In other words we have 

C(X)Q = 0; A(X) Q = a N (X) Q ; D(X) Q = 5 N (X) tt , (86) 

so that Q is an eigenstate of A(X) and D(X) simultaneously and also that for F = 
A + D. 

Other eigenvectors will be looked for in the form 

$({A}) = B(X 1 ) . . . B(Xi) ft . (87) 

The condition that $({A}) is an eigenvector of F(X) will lead to a set of algebraic 
relations on parameters Ai, ... A/. I proceed to derive these equations. 
Using the exchange relations (|67|) we get 

A(X) B(X 1 ) . . . B(Xt) Q = n /(A - A fc ) a N (X) B(X 1 ) . . . B(A,) Q + 

k=l 

I 

+ M ^ {A}) B(Ai) . . . B(X k ) . . . B(\i) B(X) Q . (88) 
k=i 



The first term in RHS has a desirable form and is obtained using only the first term 
in the RHS of (|67|). All other terms are combinations of 2 l — 1 terms, which one 
has taking A(X) to Q using the exchange relation (|67|). The coefficients M k can be 
quite involved. However the coefficient M\ is simple enough, to get it one must use 
the second term in (|67|) during the interchange of ^4(A) and -B(Ai) and in all other 
exchanges use only the first term in RHS of (^). Thus we get 

M x (A, {A}) = g(X - Ax) J[ /(Ax - A fc ) a N (X x ) . (89) 

k=2 

Now we comment, that due to the commutativity of -B(A) all other coefficients 
Mj(A, {A}) are obtained from Mi(X, {A}) by a simple substitution Ai — > Xj so that 

Mj(X, {A}) = g(X - A,) n /(A; - h) ^(X,) . (90) 

Mi 

This means, of course, that the coefficients a(A), b(X), c(A), entering i?-matrix, 
satisfy involved sum rules, making the FCR consistent. 
Analogously for -D(A) we have 

i 

D(X) B(Xt) . . . B(A,) n= l[h(X- X k ) 5 N (X) B(Ai) . . . B(X t ) SI + 

k=l 

+ £ ^(A, {A}) S(Ai) . . . B(X k ) . . . B(Xi) B(X) Q , (91) 
k=i 

where 

N^X, {A}) = k{X - A,) n HXj - A fc ) <^(A,) . (92) 

k+3 

Observe now, that 

<?(A - A,) = -k(X - Xj) . (93) 

This allows to cancel the unwanted terms in ( |S5| ) and (|5T| ) for the application of 
A(X) + D(X) to $({A}). We get, that 

(A(A) + D(A) ) $({A}) = A(A, {A}) $({A}) (94) 

with 

A(A, {A}) = o"(A) fl /(A - A,) + 8 N (X) {[ h(X - A,) , (95) 

3=1 3=1 

if the set of {A} satisfy the equations 

fl /(A, - A fc ) a^A,-) = fl h(Xj - X k ) 6 N (X j ) (96) 

k+j k^=j 

for j = 1, . . . ,1. Using the explicit expressions ( |5"9"D and ( ]55] ) we rewrite (^) in the 
form 

( Xj + i^ Y ' V-Afc+j 

U-V2J fcUAi-Afc-i • 1 j 



This is the main result of this section. In what follows we shall use the equations 
(|97D to investigate the N — > oo limit. 

An important observation is, that equations fl97|) mean, that the superficial poles 
in the eigenvalue A(A, {A}) actually cancel so that A is a polynomial in A of degree N 
as it should. This observation makes one think that only solutions {A} with Aj 7^ A& 
are relevant for our purpose. Indeed, equal A will lead to higher order spurious poles, 
cancelling of which requires more than I equations. And indeed we shall see below, 
that solutions with nonequal Aj are enough to give all spectrum. 

The equations ( j97l) appeared first (in a different form) in the paper of H.Bethe in 
1931, in which exactly the hamiltonian H was investigated. The algebraic derivation 
in this lecture is completely different from the original approach of Bethe, who used 
an explicit Ansatz for the eigenvectors $ in a concrete coordinate representation 
for the spin operators. The term Bethe Ansatz originates from that paper. We 
propose to call our approach "The Algebraic Bethe Ansatz" (ABA). The equations 
( P?D and vector $({A}) will be called Bethe Ansatz equations (BAE) and Bethe 
vector correspondingly. 

We finish this section by giving the explicit expressions for the eigenvalues of the 
important observables on Bethe vectors. We begin with the spin. 

Taking limit /i — > 00 in FCR ((P]) and using ( pE5D we get the following relation 

= , (98) 



T a (X) , \o a + S c 



which expresses the sl(2) invariance of the monodromy in the combined space 7i® V . 
From here we have in particular 



[S*,B] = -B , (99) 
[S+ , B] = A - D . (100) 



Now for reference state Q we have 



S + Q = 0, S 3 Q = y ft , (101) 



showing, that it is the highest weight for spin S c 
From ( |9"9"| ) and ( |101| ) we have 



Let us show, that 



S 3 HW) = (y - /) $({A}) . (102) 
,S + <1>({A}) = . (103) 



From ( p.00| ) we have 

S + HW) = £ S(Ai) . . . BCViX^CA,-) " J D(A,))5(A, +1 ) . . . B(X{)Q 
j 

= E O fc ({A}) J B(A 1 ) . . . B(X k ) . . . B(Xi)n (104) 



and repeating the procedure that was used to derive BAE we can show, that all 
coefficients Ok({X}) vanish if BAE are satisfied. 



Thus $({A}) are all highest weights. In particular it means, that the I cannot 
be too large, because the S 3 eigenvalue of the highest weight is nonnegative. More 
exactly we have an estimate 

N , , 

I < - . (105) 

We see that the cases of even and odd N are quite different. When N is even, the 
spin of all states is integer and there are si (2) invariant states, corresponding to 
I = N/2. For odd N spins are half-integer. 

Now we turn to the shift operator. For A = i/2 the second term (and many of its 
derivatives over A ) in A(A, {A}) vanishes and this eigenvalue becomes multiplicative. 
In particular 

Um\}) = i N F(i/2M{\}) = n^7§^({A}) • (106) 

Taking log here we see, that the eigenvalues of the momentum P are additive and 

P*(W) = E P(Ai) *({A}) , (107) 

3 

where 

P(A) = -In^| . (108) 
i A — ij 2 

The additivity property holds also for the energy H. Differentiating In A over A once 
and putting A = i/2 we get 

HH{\}) = £ *(Ai) *({A}) , (109) 

3 

where 

£ < a > = -Ia-TT74 • < 110 > 



Formulas (|108|) and (|110|) allow to use the quasiparticle interpretation for the 



spectrum of observables on Bethe vectors. Each quasiparticle is created by operator 
5(A), it diminishes the S 3 eigenvalue by 1 and has momentum p(X) and energy e(A) 
given in (|108|) and (|110|) . Let us note, that 

£ (A) = i^p(A). (Ill) 

The variable A in this interpretation can be called a rapidity of a quasiparticle. 

It is possible to exclude the rapidity to get the dispersion relation, describing 
connection of energy and momentum 

e(p) = cosp — 1 . (112) 

The eigenvalues of hamiltonian are all negative, so that the reference state Q 
cannot be taken as a ground state, i.e. state of the lowest energy. It trivially 
changes if we take —H as a hamiltonian. Both cases H and —H are interesting 
for the physical applications, corresponding to antiferromagnetic and ferromagnetic 
phases, correspondingly. The mathematical (and physical) features of the N — » oo 
limit in these two cases are completely different, as we shall see soon. 



5 XXXx/2 model. Physical spectrum in the fer- 
romagnetic thermodynamic limit 



In our case the thermodynamic limit is just limit N — > oo. We shall see, how BAE 



simplify in this limit. Looking at BAE ( 97|) we see that N enters there only in the 



exponent in the LHS. For real Aj, . . . , Aj both sides in BAE are functions with values 
on the circle and LHS is wildly oscillating when N is large. Taking the log we get 

i 

Np(Xj) = 2nQ j + £ Lp(\j - A fc ) , (113) 

k=l 

where the integers Qj, < Qj < N — 1 define the branch of the log and <f(X) is 
a fixed branch of ln^r. For large N and Q and fixed I the second term in the 
RHS of ( |113j ) is negligible and we get the usual quasicontinuous expression for the 
momentum of a free particle on a chain 

Pj = 2vr^ . (114) 

In the ferromagnetic case when the hamiltonian is —H the energy of this particle is 
given by e(p) = 1 — cos p. 

The correction (second term in RHS of ( |113j )) expresses the scattering of these 
particles. The comparison with the usual quantum mechanical treatment of a parti- 
cle in a box shows that <p (Aj — A&) plays the role of the phase shift of particles with 
rapidities Aj and A*.. Thus the function 

S(A-Ai) = ^±i (115) 
A — p, — % 

is a corresponding S'-matrix element. 

Another analogy is with the classical inverse scattering method, where the com- 
bination 

Z(X) = 5(A)A- X (A) (116) 
was more important, than B(\). The factor S enters the exchange algebra 

Z(\)Z(ji) = Z(ji)Z(\)S(\-fi), (117) 



valid in the limit iV — > oo because the second term in the RHS of (p7\ ) effectively 
vanishes. Operator Z(\) can be interpreted as a creation operator of a normalized 
particle state. 

This argument by analogy should be justified by the adequate scattering theory 
applicable to our case. We do not have time to do it here and refer to the original 
papers of Babbit and Thomas. 

The BAE allow also for the complex solutions which in our situation correspond 
to bound states. Let us see it in more detail. The first nontrivial case is I — 2. From 
two BAE 

{X^ij2 " Al -A 2 -i ' (118) 



(120) 



( X 2 + i/2 \ N A 2 -A 1 + z 

we see, that 

/ A 1 + z/2 \ jV / A 2 + z/2 \ Af _ 
\Xi - i/2 J \X 2 -i/2) 

so that p(Ai) + p(A 2 ) is real. Further, for ImAi 7^ the LHS in ( |118| ) grows (or 
decreases) exponentially when N 00 and to compensate it in the RHS we must 
have, that 

Im(Ai - A 2 ) =i (or - i) . (121) 

As Ai and A 2 can be interchanged, we can say, that in the limit N — > 00 Ai and A 2 
acquire the form 

t % 
Ai = A1/2 + - , A 2 = A1/2 - - , (122) 

where Ai/ 2 is real. In the thermodynamic limit Ai/ 2 can become arbitrary. 

The momentum pi/ 2 (A) and energy ei/ 2 (A) for the corresponding Bethe vector 
are given by 

JP1/2W = j P0 (x+i/2)+ip (x~i/2) = X + 1/2 + i/2 _ A + i/2 - t/2 = \ + i 

A - i/2 + i/2 X- i/2 -i/2 X-i y J 

and 

£l/2W = ii_, npi/2(A) = _L_. (124) 

The origin of notation po(A) for the momentum ( |108| ) and Pi/ 2 (A), ei/ 2 (A) for mo- 
mentum and energy of the complex solution will be clear soon. Excluding A from 
(|123D and (ggg) we get 

e 1/2 (p) = i(l-cosp). (125) 
The interpretation of this eigenvector as bound state is supported by the inequality 

ei/ 2 (p) < e (p-pi) +eo(pi) (126) 

for all p, pi, < p, pi < 27T. The S'-matrix elements for scattering will be discussed 
later. 

For I > 2 the complex solutions are described analogously. Roots A; are com- 
bined in the complexes of type M, where M runs through half-integer values 
M = 0,l/2,1,... , defining the partition 

Z = $>m(2M+1) , (127) 

AI 

where Vm gives the number of complexes of type M. 

The set of integers {vm} defines a configuration of Bethe roots. Each complex 
contains roots of the type 

X M , m = A M + im , -M < m < M , (128) 



where Am is real, m being integer or half-integer together with M. The momentum 
and energy of complex is given by 



1 A + z(M + l/2) . , 

and 

/A , 1 2M+1 1 , . . 

£m(a) = = ^TT (1 " c ° SPm) " (130) 

The S-matrix element for the scattering of complex of type on a complex of type 
M is given by 

\ + iM \ + i(M + l) 
SoMX) ~ A^M ' A^mTI) ( } 

and for scattering of complexes M and N 

M+N 

Sm,nW= II 5 o,l(A). (132) 

i=|M-JV| 

It is the superficial analogy of this formula with the Klebsch-Gordan formula for 
sl(2) which prompted me to use label M for complexes instead of their length 2M+1. 
The derivations are just the direct calculation of products 

, ' X + i/2 + im X + i + im ^ A + i + i{m + n) 

where many terms cancel. 

This finishes the description of the physical spectrum of —H in the thermo- 
dynamic limit. The ground state Q = Y\u n defines the incomplete infinite tensor 
product in the sense of John von Neumann. The space Tip is a completion of the 
states which differ from u n only in finite number of factors in Y\®h n . The excita- 
tions are particles, classified by half-integers M, M = 0, 1/2, 1, . . . and rapidity A 
(or momentum p) . The dispersion law for a particle of type M is given by (|130|) and 
scattering matrix elements are given by (|132|). The interpretation of particles with 



M > as bound states is possible but not necessary. 

The spin components S have no sense in the physical Hilbert space 7Yf as they 
change vector Q in any site n. Operator S 3 after shift by N/2 

Q = f - S 3 (133) 

makes sense in Tip and has integer eigenvalues 2M + 1 . This phenomenon gives the 
example of "symmetry breaking" by vacuum: only U(l) part of sl(2) remains intact 
in the thermodynamic limit. Physically we call this phase ferromagnetic because 
the prescribed direction of spin (in 3-d direction in spin space) is macroscopic and 
fixed. 



6 XXX x / 2 model. BAE for an arbitrary configu- 
ration 



Before turning to physics of the antiferromagnetic chain we shall consider in more 
detail the BAE for arbitrary configuration {vm}, where each integer vm gives the 
number of complexes of type M combined from / quasiparticles, so that 

Z = + l)i/ M . (134) 

M 

We shall investigate the approximate BAE for the real centers of complexes Am,;, 
and allow I to be of order N/2. There are some doubts, expressed in the literature, 
about validity of the picture of complexes in this case. Indeed, in our arguing above 
we supposed, that I is much smaller than N. However in the physical applications 
only uq will be large and in this case the picture of complexes is correct. 

The BAE for the Xmj, j = 1, • • • , are obtained by multiplying the BAE for 
each complex in the LHS of (p7|) and rearranging the RHS according to the picture 
of complexes. They look as follows 

e ipM(A ^ )7V = II II S M , M ,{\ Md -\ M ',k). (135) 

M' (M',k)^(MJ) 



The factors in the RHS are scattering matrix elements from ( |132| ); in the LHS the 
momentum ^^(A) from ( |129| ) enters; condition (M',k) ^ {M,j) means, that among 
vm> roots \M',k, entering the RHS in ( |135| ), the one which is equal to Xmj from LHS 
is absent. 

Taking the logarithm of (|135|) and using the basic branch in the form 



- hi ^ =7r — 2arctg— (136) 
i A — ia a 



we get the equation 

Xmj 



2ATarctg * = 2nQ Md + ]T £ <I\y..w<( A.w., - \ M >,k) , (137) 



where 



M+M 1 



$m,m'(A) = 2 ( arctgy + arctg-^— j (138) 

L=\M-M'\\ L 

with the understanding that the term with L = is omitted and Q mj is an integer 
or half-integer (depending on the configuration, which parametrizes the roots Xmj)- 

The main hypothesis of our investigation is that Qmj classify the roots uniquely 
and monotonously: roots Xmj increase, when Qmj increase; moreover there are no 
coinciding Qmj for a given complex of type M. The last condition corresponds to 
the requirement for the BAE roots to be distinct, which was discussed in the course 
of derivation of BAE. 

We shall look for the real and bounded solutions of equation ( |1371) . For those 
the numbers Qmj have a natural bound. Indeed, taking into account, that 

arctg±oo = ±^ (139) 



and putting Xmj = oo we get for the corresponding Qmj the expression 

Qm,oo = - £ (2 min(M, M') + 1) - (W + J) - 1) + ^ . (140) 



The maximal admissible Qm,j is then 

QM,max = Q A f,oo ~ (2M + 1) (141) 

because complex of type M has (2M + 1) roots. We suppose, that when Qmj gets 
values bigger than QA/.max, the roots in our complex turn to be infinite one after 
another, so that for Qmj = Qm,oo the whole complex becomes infinite. 

I understand, that all this is quite a host of hypotheses, but the result we shall 
get soon is quite satisfactory. It will be nice to produce more detailed justification 
for our considerations. 

From ( |140| ) and ([HI]) we get 

QAf,max = ^-E J{M,M')v M <- \ , (142) 
1 M> 1 



where 



T(M m>\-\ 2min(M,M') + l M + M> 



2M+| M = M' . 



(143) 



Analogously we find QM,mm- Due to the fact, that arctgA is odd we have 

Qjtf,min = —QM,max ■ (144) 

Thus for the number of vacances Pm for the numbers Qmj we have 

Pm = 2g A f,max + 1 = N - 2 J(M, M')v M , ■ (145) 

w 

The numbers Qmj are integers for odd Pm and half-integers for Pm even. 

Now we can estimate the number of Bethe vectors, characterized by the admis- 
sible numbers Qmj- For a given configuration {vm} the states are given by fixing 
the distribution of Q-s over the vacances; so the whole number Z(N, {vm}) of them 
is given by 

Z{N } {u M })=\{CZ , (146) 

M 

where is a binomial coefficient 

in ' 

m\[n — my. 

Let us consider the number of states for given I and number of complexes 

q = Y, v M (148) 

inside each configuration {vm} 

Z(N,l,q)= Z{N-{v M }) . (149) 



We shall calculate Z(N, I, q) by reduction via a partial summation. For that we 
shall begin by extracting the contribution of roots of type 0, or in other words by 
substituting the configuration {vm} by {^m}, where 

z4r = ^M+i/ 2 , M = 0,^,... . (150) 

First we observe, that 

Pm(N, {u m }) = P M _ 1/2 (N - 2q, {u' M }) . (151) 
Indeed, it is easy to see that 

J(M, M') = J(M - i M' - ~) + 1 , (152) 

so that for M > 1 /2 

P A f(iV,{z/ M }) = iV-2J(M,0)z/o-2 ^ J(M,M')v M ,= 

M'>l/2 

= N-2u -2 ]T (j(m-~, M'-^+l\v M > = 



N-2q-2Y / j(M-l,M')u' M , (153) 



M'>l/2 

/f -2' 

A/' x Z 

and ( |151| ) follows. Thus we have a recurrence relation 

Z(N, {u M }) = C%Z (N - 2q, {v' M }) (154) 
and summing over the allowed vq we get 

Z(N, l,q) = J2 C^ 2q+U Z(N -2q,l-q,q-u). (155) 
With the initial condition 

Z(N, 1, 1) = JV - 1 (156) 

this gives 

W ?) = ^" 2 /^ 1 1 C^_ i+1 C?_ 1 (157) 
and finally for the number of the Bethe vectors with given / 

Z(N, I) = £ Z, g) = C l N - C^ 1 . (158) 

9=1 

Now we remember that each Bethe vector of spin y — Z is a highest weight in 
the multiplet of dimension N — 21 + 1. Thus the full number of states, described in 
our picture 

Z = J2(N - 21 + 1)Z{N, Z) = 2 N (159) 
i 

is equal to the dimension of our Hilbert space. This is very satisfactory and strongly 
confirms all the hypotheses, which we used in this calculation. We stop here the 
general investigation of the BAE (|137|) . 



7 XXX x / 2 model. Physical spectrum in the anti- 
ferromagnetic case 

I am ready now to describe some important states. The ground state, i.e. the 
state of the lowest energy, is obtained by taking the maximal number of real roots. 
We shall suppose, that N is even to be able to have an si (2) invariant state. The 
corresponding configuration looks like 

N 1 
^o = y; ^m = 0, M>-. (160) 

For this configuration I = N/2 and so 

S 3 = --l = 0; (161) 

thus the spin of the state vanishes. Now the number of vacances P 

N N 

P = N- 2.7(0, 0)v = N- — = — (162) 

is equal to the number of roots and so there is no freedom for the allocating the 
numbers Qo t k] they span all interval 

N 1 N 1 

- j + ^ < Qo,t < j - ^ (163) 

being integer (half-integer) for N/2 odd (even). 

Thus the state in question is unique and gives us the singlet for sl(2) group of 
spin observables. 

Using physical terminology we can call this state the Dirac sea of quasiparticles. 
The mere existence of this state is due to the Fermi character of the quasiparticles 
spectrum, i.e. to the condition that all Qqj are distinct. 

In what follows we consider states, for which u differs from its maximal value 
N/2 by a finite amount 

N 

u = j-k . (164) 

We shall see, how a Fock-like space of excitations will emerge step by step with 
increasing of k in the limit iV — > oo. Thus it is k, which will play the role of 
the "grading" in our definition of the physical portion in the formal infinite tensor 
product n®C 2 . 

For a fixed k all v M , M > 1/2 are bounded, when iV — > oo. Indeed, from 
inequality I < N/2 it follows 

£ (2M + \)v M = 1-v <^-v = k (165) 

M>l/2 



and N disappears from this estimate. Thus for a fixed k the number of configurations 
is finite and we can consider them one after another. 



For k — 1 only um = 0, M > 1/2 are allowed. The / for these states is I = N/2—1, 
thus the spin is 1. The number of vacances 



1 /N \ N 

P = JV _ 2 ._(__ 1 ) = _ + 1 (166) 

exceeds the number of roots by two. Thus two admissible numbers Qqj are not to 
be used, which gives the two parameter degeneracy of the state. In physical jargon 
one speaks of the holes in the Dirac sea. 

For k = 2 there are two possibilities: vm — 0, M > 1/2 and v\/2 = 1, vm = 0, 
M > 1. Let us consider the latter in more detail. First, I = N/2 for it, so that the 
spin vanishes. Second, for the corresponding Pq we have 

fl, = tf-2(£-2).I-2j(0,I)=£ (167) 

and 

^ = W -2(f-2)j(i,0)-2j(i,i) =4-3 = 1. (168) 

We see, that the number of vacances for real roots once more exceeds their number 
by 2 and there is no freedom at all for the root of complex of type 1/2. 
In the former case 

P --_2(^-2)I^ + 2 (16, 

exceeds the number of roots by 4 and the spin of the state is 2. 

For general k we have a similar picture. First, for the configurations 

N 

u = --k; v M = 0, M > 1/2 (170) 

we have 

N 

Po = y + k > ( m ) 

so that there are 2k holes, characterized by the missing places in the choice of 
admissible Qo,j- The spin of this state is equal to k. Then there are states with a 
smaller spin, corresponding to a few nonzero i>m, M > 1/2. More on this will be 
said later. 

We return to a more detailed characteristic of the states already described. For 
this more control over the roots of BAE is needed. Fortunately these equations are 
simplified drastically in the thermodynamic limit iV — > oo. The real roots become 
quasicontinuous and we can evaluate their distribution. 

We begin with the ground state. The roots are real and corresponding Q j 
fill without holes all interval (|163|) . We shall put for N/2 odd (with an evident 



correction for N/2 even) 

Qo,j = J , (172) 

so that the BAE take the form 

arctg2Aj = j- + — ^ arctg(A j - \ k ). (173) 

k 



The variable 

x = jj (174) 

becomes continuous in the limit iV — > oo with values —1/4 < x < 1/4; the set of 
roots Xj turn into function X(x). 



The equation (|173|) becomes 

,1/4 

arctg2A(x) = ttx + / arctg (A(x) — X(y)) dy (175) 

7-1/4 



and looks rather formidable. Fortunately it is not X(x) which is of prime concern to 
us. Indeed, we are interested in the eigenvalues of local observables, which take the 
form of the sums over roots (see e.g. ( 107[ )). With our conventions we have 



/1/4 roo 
h(X(x))dx = N / h(X)p(X)dX , (176) 
-1/4 J-oo 



where the change of variables A : x — > X(x) maps interval — 1/4 < x < 1/4 into 
whole real line — oo < A < oo due to the monotonicity of X(x). The density p(A) is 
nothing but 

dx 1 
p( X ) = IT 



dX X'(x) 



(177) 

x=X- 1 (X) 



Differentiating (|175|) we get for this density, which we denote by po(X) for our 
state, a linear integral equation 

WA) + / ~j ' 7t ^dji , (178) 



1 + 4A 2 ruv ' J-oo l + (A-/i) 2 
which can be easily solved by Fourier transform. We get 

The momentum and energy of the ground state are given by 

Po = N J Po (X)p (X)dX = (180) 

due to the fact, that the integrand is odd if we use for p (A) slightly shifted expression 

p (X) = -2arctg2A (181) 

and 

E = N e (A)p (A)dA = -iVln2 . (182) 

Note, that the sign of e^,/(A) in this section is opposite to that in section ^. Thus 
the energy of the ground state is proportional to the volume as always in the correct 
thermodynamic limit. Adding it to the hamiltonian H will make it nonnegative. 

Now we turn to the configuration v$ = N/2 — 1, vm = 0, M > 1/2. There are 
two holes and we can put 



Qoj =3 + 0(3 ~ h) + 0(3 ~ 32) 



(183) 



where 6 is a step function 



e<j) = \l HI d84) 



and ji and j 2 are integer, characterizing the holes. By the same trick as before we 
get for the distribution pt(A) (t for triplet) of real roots the linear integral equations 

^t(A) + [°° - , P A (/i) ^ dfj, + ^(5(X - Ai) + 5(X - A 2 )) , (185) 



1 + 4A 2 r yj J -00 1 + (A - /i) 2 N 

where Ai and A2 are images of x\ = ji/N and xi = j 2 /iV in the map x — > X(x) 
defined by A t (x) (or A (x), because Ai and A 2 enter in terms of order l/N). From 



185| ) we get 

Pt(A) = p (A) + ^(a(A - Ax) + a{\ - A 2 )) , (186) 
where cr(A) solves the equation 

ff(A) + 7/lTTxr^* + * (A) = °- (187) 

Solving it we can evaluate the momentum and energy of the corresponding state 

P =k(\i) + k{\ 2 ) ■ (188) 
E = E + fc(Ai) + MA2) , (189) 

where 

MX) = arctgsinhvrA, e(A) = ^ — -. (190) 

2 cosli7rA 

It is time to comment, that the constructed states once more allow for the particle 
interpretation: we have described a family of two particle states with the energy and 
momentum of a particle given by ( |190| ) and dispersion law is 

e(jfe) = |cosA;, -tt/2 < k < tt/2 . (191) 

Next example is u Q = N/2 — 2, vi/2 = 1, vm = 0, M > 1. For the density of real 
roots p s (A) (s for singlet) we get the equation 



Tip, 



(A) + looiTTA^^ + 



1 + 4A 2 rsv ' J -00 l + (A-p) 

+ 1 (5(A - Ai) + 5{X - A 2 ) + $' 0il/2 (A - A 1/2 )) , (192) 

where Ai and A 2 stand for the holes and the last term in the RHS is a contribution 
of the complex of type 1/2. For Ai/ 2 we have one more equation 

arctgAi /2 = — $i/2,o(Ai/ 2 - X 0J ) , (193) 
i 

because the corresponding number Q1/2J has just one admissible value equal to zero. 



From ( |192| ) we get 



p s (A) = p (A) + -(a(A - Ai) + a(X - A 2 ) + cu(A - A 1/2 )) , (194) 
where cr(A) is as above and u is a solution of equation 

/oo Cjf/i") 
_ i + ( A-^ rf/i + *M A ) = • ( 195 ) 

To evaluate A 1/2 let us rewrite ( |193| ) in the form 

/oo 
$i/2,o(Ai/ 2 - A)p (A)ciA = (196) 
-oo 

1 f 00 

= -77 / $i/2,o(Ai/ 2 - A) [cr(A - Ai) + <r(A - A 2 ) + w(A - Ai /2 )]dA , 

iV J— 00 

where the limit iV — > 00 is already taken into account by changing the sums over 
Aoj by integral over density p B (X). LHS here vanishes for any A1/2 and contribution 
of w disappears due to oddness of integrand ; so we get the equation 

/oo 
$i/2,o(Ai/ 2 - A) (a (A - A x ) + <r(A - \ 2 ))d\ = , (197) 
-00 

or 

arctg2(Ai /2 - Ai) + arctg2(Ai /2 - A 2 ) = (198) 

with the solution 

A 1/2 = ^ • (199) 

We are ready now to evaluate the observables. The spin of our state is zero. For 
momentum and energy we get exactly the same expressions ( |188| ) and ( |189| ) as in 
the previous example; the contribution of a complex of type 1/2 cancels exactly. 

The examples considered are all, which give a two-parameter family of states. 
Returning to the particle interpretation, we can say, that our particles have spin 
1/2. Indeed, we constructed the highest weights in triplet and singlet two-particle 
states with exactly the same momentum and energy content. Thus they are highest 
weights in C 2 <S> C 2 representation of a spin observable. This is why I say, that the 
particles have spin 1/2. 

This picture is recurrently confirmed in the description of the next excitation. 
The state is = N/2 — k, I'm — 0, M > 1/2 defines a 2k particle state being the 
highest weight in the highest spin irreducible component in Y[ 2k ®C 2 . All other states 
for the same k are states of lower spin, entering into multiplets with the number of 
particles not exceeding 2k. The contribution of complexes of type M into energy 
and momentum always vanishes, so that the energy-momentum expressions depend 
only on the number of particles. 

All this allows to say, that the only excitation of our system is a particle with spin 
1/2 and energy - momentum relation ( |191| ). Note, that the momentum k(X) runs 
through the half of the usual Brillouin zone. There is one important restriction: the 
number of particles is even. These particles are usually referred to as spin waves. 
For a long time it was stated in the physical literature, that spin waves of the 



antiferromagnetic chain of spin 1/2 magnets has spin 1. Indeed, spin wave being a 
hole in the singlet Dirac sea corresponds to a turn of one spin, amounting to spin 
1/2 + 1/2 = 1. 

However, our more precise analysis shows, that turn of a spin corresponds to 2 
holes and two spin wave excitations. The momentum of this state runs through the 
whole Brillouin zone —it < k < tt. 

Mathematically the Hilbert space Haf is a (half) Fock space 



•T/CVCIl 
ft 



00 rir/2 ptt/2 2n 

AF=E/ dm... dn 2n \{®C 2 . (200) 

J-tt/2 J-tt/2 



n=0' 



t/2 J-tt/2 



(Compare it with Professor's Miwa lectures, where TCaf corresponds to A <8> A , or 
Ai <g> Ai). 

Natural question is how to describe one particle (or odd number of particles) 
state. The answer is, that they enter the chain of odd length. The lowest energy 
state there has spin 1/2 and have one hole in the distribution of numbers Q ,j- 
Thus this state becomes a one-particle state in the thermodynamic limit. Chain 
of odd length has no ground state but just 1 particle, 3 particles etc. states. The 
corresponding Hilbert space is 



oo r n/2 2n + 1 

HIV = E / d Kl ... dK 2n+l J] ®C 2 (201) 

n=0 J - 7T / 2 J-tt/2 



(and corresponds to A <8> Ai or Ai <8> A in Professor's Miwa lectures). 

We finish with some formalization of our result. The expression for the Bethe 

state 

mX}) = B(X 1 )...B(\ l )n (202) 

can be formally rewritten as 

$ ({A}) = jexp £ lnBM j SI (203) 

and now it is possible to go to the thermodynamic limit. For the ground state $ 
we have 

$o = exp|^ J lnB(\)po(\)d\^n . (204) 

Exponential dependence of a true ground state on the volume is a typical phe- 
nomenon in quantum field theory. Now the triplet excited state $(Ai,A 2 ) can be 
written as 

$(Ai,A 2 )=Z(Ai)£(A 2 )$ (205) 
without any reference to Q and dependence on the volume. Here 

Z{\) = exp y In B(X)a(X - /i)rf/i| (206) 



plays the role of a creation operator of one-particle state from the physical ground 
state. 



As was mentioned in the section ^ it is 

Z{\) = Z(\)A-\\) (207) 

which is a more natural object in the scattering problem. Of course Z(X) is defined 
up to a constant normalization factor, which we can use to cancel the iV-dependent 
factor cioo(A) 



«oo(A) = (A + ^) Jv exp|iV J_Jn \^ * PofaWj (208) 
entering the eigenvalue 

A N (\)Q = aoo(A)$ (209) 

in the limit N — » oo, when the contribution of D^(X) to A(A, {A}) vanishes. 
Operators Z(\) satisfy the exchange relation 

Z(X)Z(jj,) = Z(n)Z(X)S t (X - fj,), (210) 

where the phase-factor S t (A — /x) is given by 



St (A) = exp <( / In- ; <7(/i- A)d/i^ = y ^Ti^XVFTi , flr (2II) 



oc 



In course of derivation the contribution of the second term in the exchange relations 
(|67D is negleqted, which can be justified in the limit iV — ► oo. 

The factor S t (A) is to be interpreted as a triplet eigenvalue of the S-matrix for 
spin 1/2 particles, acting in C 2 ® C 2 

S^\X-ri = S t (X-ri(-^-I+-^—p) . (212) 

\X — n + i X — fi + i J 

The creation operators Z e (X), e = ±1 are to satisfy the exchange Zamolodchikov 
relation 

Z ei {X x )Z e2 {X 2 ) = Z 4 (A 2 )Z 4 (A 1 )<^(A - V) ■ (213) 

We did not construct operator Z £ (X) (the vertex operators of the second kind in 
Professor Miwa terminology). We can only identify in 7i c ^§ n 

Z + (X)Z + (») = Z(X)Z(») , (214) 

where Z(X) is given by (p07|) and 

Z+(\)Z-(») - Z-(\)Z + {») = (215) 



OO 



Z(X) exp I I In B(a)u ( *±± - a ) da] B ( A + ? + M 5 ( A+ ?"M Z^). 



— oo 



This identification is, however, sufficient to justify all S-matrix 

The interesting but not understood comment on the formula ( [21 1| ) is as follows: 
the phase-factor St (A) coincides with the S-matrix for the rotationally symmetric 
subspace of the Laplacian on Poincare plane. Indeed putting s = (1 + iX)/2 we get 



where 

= roT^TT) (217) 

is a Harrish-Chandra factor for sZ(2, R). With this intriguing comment we finish 
our long and detailed treatment of the XXX1/2 model. From now on I shall describe 
several directions of development and/or generalization along the similar lines with- 
out giving too much details. The first generalization is a XXX model for higher 
spin. 



8 XXX S model 

Now I consider the spin chain with local spin variables realizing the finite di- 
mensional representation of si (2) in 2s + 1 dimensional space C 2s+1 , where s is spin, 
integer or half-integer. I am not ready to write the corresponding hamiltonian. To 
maintain the integrability I shall find it as a member of the commuting family of 
operators, generating function for which will be trace of an appropriate monodromy 
of the family of local Lax operators, satisfying the FCR a-ld ( 44]) . 

The Lax operator L n a (A) with the auxilialy space V = C 2 does not differ from 
(j32|) . Indeed, operator, defined in h n ® V = C 2s+1 <g> C 2 by matrix 

L n , a (X) = XI + • £ Sy a = ( A + f » A ^ ) (218) 

satisfy the relation 

(A-//) (219) 

with the same i?-matrix R aia2 (X) from (^Sj). The derivation from section |3] is not 
applicable. We shall not derive ( |219|) here because a more general check will be done 
below for the XXZ model. 

Introducing the monodromy 



T a (X) - II^n,a(A) - ( 0^ j 



(220) 

we see, that it satisfies FCR of the form (f44|) , so that its trace 

F(X) = A(X) + D(X) (221) 

is a commuting family of operators 

[F(A),F(/i)] = 0. (222) 

This family can be diagonalized by means of Algebraic Bethe Ansatz (ABA). Indeed, 
we have local vacuum uj n - the highest weight in C 2s+1 , the reference state 

Q = Y[ ®uj n (223) 

with the eigenvalues for A(X), -D(A) and C(A) 

A(X)n = a N {X)Q, D(X)n = 5 N (X)n, C{X)Q = , (224) 



where 

a(X) = X + is, 5(A) = X -is, (225) 

and exactly the same exchange relation for A(X), -D(A) and -B(A) as (|67|)-(|68|). Thus 
the state 

m\}) = B(\ 1 )...B(\ l )tt (226) 
is an eigenstate of the family F(X) with the eigenvalue 

A(A, {A}) = (A + is) N {[ + (A - is) N {[ , (227) 

j= i A - A j i= i A - A j 

if {A} are roots of the BAE 

N 



fXk±is\ =n V-A I±1 

VA fc -^y ^Xk-Xj-i 

When s = 1/2 we return to the case already considered above. 

However the construction of the local hamiltonian cannot repeat one from section 
||. Indeed, in no point A the Lax operator L n)Q (A) is a permutation, just because the 
quantum space h n = C 2s+1 and auxiliary space V = C 2 are essentially different. 

The way out is to find another Lax operator, for which the auxiliary space V is 

Q2s+1 

The existence of such an operator is based on a more general interpretation of 
the FCR due to V.Drinfeld. 

In this interpretation the generating object for Lax operators is a universal R- 
matrix 1Z defined as an element in A® A for some algebra A, satisfying the abstract 
Yang-Baxter relation (YBR) 

n 12 n 13 n 23 = n 23 n 13 n 12 . (229) 

This equation holds in A ® A <E> A and rather evident notations are used, i.e. 

7^12 = n ® 1, n 23 = \®n . (230) 

The algebra A must have a family of representations p(X,a), parametrized by a 
discrete label a and continuous parameter A. For instance the loop algebra of any 
unitary group has such representations called the evaluation representations. In our 
case of XXX models the algebra A was identified by Drinfeld and called Yangian 
by him. Below on the case of XXZ models we will encounter the g-deformed affine 
algebra as A. 

The concrete Lax operators are obtained via the evaluation representations of 
the universal 7£-matrix, i.e. 

L n , a (X-fi) = (p(a,X)®p(n,fi))1l = R a ' n (X-p) . (231) 



The dependence in the LHS on A — p reflects some homogeneity in the family of 
representations p(a, A). The Yangian for si (2) has representations p(a, A), where a 
is just spin label of representations of sl(2), a = 0, 1/2, 1, . . . . The relation (|36D is 



obtained if we apply the representation p(l/2, A) ®p(l/2, /i) ®p(s, cr) to YBR ( |229| ), 
put <T = and identify 

# ai>a2 (A) = R 1/2 ' 1/2 W ; ^n,a X (A) = R'^iX). (232) 

However we can use another combination of the representations. Let us rewrite the 
YBR (|229|) in the form 

n l2 n 32 n 3l = n 31 n 32 n 12 , (233) 



which can be easily derived from ( |229| ) by applying the appropriate permutation to 
it. Now apply to (|233|) the representation p(s\, A) (g) p(s 2 , fi) ® p(l/2, cr). We get 

R S ^ S2 (X - p)R 1/2 ^(a - p)R 1/2 ' Sl (a - A) = 
= R l ' 2 ^{a - X)R 1/2 ' S2 {a - (i)R Sl ' S2 {X- /i) . (234) 

The factors i? 1//2,s (A) can be identified with the Lax operators L n a (A) above; the 
operator R Sl ' S2 (X) give us a new Lax operator; we can take representation with spin 
s± as a local quantum space and that with spin s 2 as an auxiliary space. In particular 
for s\ = s 2 we get the Lax operator L n j(\) we are looking for. Equation ( [234j ) is a 
linear equation, which allows to calculate L n j if L n>a are known. We shall call L n j 
the fundamental Lax operator and label / stems from this. 

Let us now calculate R Sl ' S2 (\) for representations s x and s 2 being the same, using 
equation ( |234[ ) . To simplify the notation we shall denote two sets of spin variables 
by S a and T a and use the notations 

L T (X) = \ + i{T,a), L S (X) = \ + i(S,a) (235) 

for the corresponding Lax operators i? 1 / 2 ' 5 . Here 

(T,a) = J2T a a a (236) 

a 

and analogously for (S, a). We shall look for R Sl,S2 (X) in the form 

R S ^ S2 (X) = P Sl ' S2 r({S,T),X) , (237) 
where P S1,S2 is a permutation in C 2s+1 ® c 2s+1 and {S,T) is a Casimir C 

C = {S,T)=Y, S a T a . (238) 

a. 

Using 

P(S»P = (7» (239) 
we rewrite the equation ( |234| ) as follows 

(A - i(T, <r))(ji - i(S, a))r{X - fx) = r(A -fi)Qi- i(T, a))(A - a)) . (240) 

We have due to the property of Pauli matrices a a 

(T,a)(S,a) = (T,S) + t((SxT),a) , (241) 



where 

(SxT) a = e a ^r . (242) 
Now using the central property of Casimir we transform equation ( |240f) into 

{XS a + (T x S) a )r{\) = r{X){XT a + (T x S) a ) . (243) 

Due to symmetry, it is enough to consider one out of three equations ( |243| ) e.g. the 
combination 

{XS + + i{T 3 S + - S 3 T + ))r{\) = r(A) (AT+ + i{T 3 S + - S 3 T + )) . (244) 

We shall use a convenient variable J instead of Casimir (S,T); taking into account 
that representations for 5* and T are irreducible we have 

(S + T) 2 = S 2 + T 2 + 2{S, T) = 2s(s + 1) + 2{S, T) = J{J + 1) , (245) 

where the operator J have an eigenvalue j in each irreducible representation Dj 
entering the Klebsch-Gordan decomposition 

2s 

D s ®D s = Y, D o ■ ( 246 ) 

3=0 

We shall look for operator r(A) as a function of J. To find it we shall use the 
equation ( [244j ) in the subspace of the highest weights in each Dj, i.e. put 

T + + S + = . (247) 

This is permissible because 

\T + S 3 - S + T 3 , T + + S + ] = . (248) 
In this subspace due to the fact, that in general 

(S + Tf = {S 3 + T 3 ) 2 + S 3 + T 3 + {S- + T~)(S + + T+) (249) 
we can identify 

J = S 3 + T 3 . (250) 
In the constrainted subspace the equation (|244|) reduces to 

(XS + + tJS + )r{\, J) = r(A, J)(-XS + + tJS + ) (251) 

Now we use the commutation relation 

S+J = S + (S 3 + T 3 ) = {S 3 + T 3 - 1)S + = ( J - 1)S + (252) 



to turn (|25 1|) into the functional equation 

(A + iJ)r(\, J - 1) = r(A, J)(-A + U) (253) 

with solution 



normalized in such a way 

r(J,0) 



r(J,-\)r(J,\) = l . 



(255) 



Of course in our case of finite dimensional D s we are interested only in r( J, A) for J 
taking values J = 0, 1, 2, . . . , 2s. But in what follows we shall use the representations 
with any complex s and then formula (|254j) will be used in all generality. 
Having this Lax operator L n j(X) we can write one more variant of FCR 

R h,h(^ ~ ^LnjtWLnj^) = L n j 2 ([i)L nJl (\)R fl j 2 (\ - /i) . (256) 

From this we infer, that the spectral invariants of the monodromy 

T f (X) = L NJ (X)L N ^j(X) . . . Lij(A) (257) 

(i.e. Ff(X) = tifTf(X)) are commuting 

[F f (\),F f fa)] = . (258) 

The relation ( |234| ) with si = S2, where i? 1/,2,s (A — /i) is used as i?-matrix and i? 1//2 ' s (^) 
and R S ' S (X) as Lax operators leads to the commutativity of the families Ff(X) and 

[F f (\),FM] = . (259) 

Thus we can use Ff(X) to get observables and F a (X) to construct BAE. 
Repeating the considerations in section ^| we get 



F f (0) = U = e 



iP 



and 



where 



H 



d 



N 

= E Hn, 

A=0 n=l 



71+1 



H n> n+x = i — lnr(J, A) 



(260) 
(261) 
(262) 



A=0 



where J is constructed via local spins and S^_^^ as 

J( J + 1) = 2 J2(S:S: +1 ) + 2s(s + 1) (263) 

a 

From (g5|) and (|262|) we get 

Hn,n+1 = —2ijj{J + 1) , (264) 
where ip(z) is logarithmic derivative of T(z). For positive integer n we have 

# + n ) = El-7 > (265) 
fe=i K 

where 7 is the Euler constant, and it allows to express H n ^ n+ i as a polinomial in 
invariant J2 a ^n^n+i = C n ,n+i as follows 



-f^n,n+l — E Cfc (^ n .™+ 1 ) fc ~~ f2s(C n ,n+l) 



(266) 



where the polinomial f2s{x) can be written using Lagrange interpolation as 



/*(*) - E ( E \ - 7 ) n fnv' x < = ^ (/(/ + 1} - 2s(s + 1)} • (267) 

In particular for s = 1 we have 

ci = -c 2 , (268) 

so that the hamiltonian 

H = £ - (S:^ +1 ) 2 ) (269) 

is integrable for the representation of local spins in C 3 . The naive generalization of 
the hamiltonian ( p7|) by simple substitution of operators of spins 1/2 by those of 
spin 1 is not integrable. 

The construction of the integrable hamiltonians for spin s magnetic chains is one 
of real achievements of the Algebraic Bethe Ansatz. Indeed, without well under- 
stood connection of integrable hamiltonians and Lax operators there is no hope to 
reproduce the formulas ( |266|) , ( [26 7j ). 

Now I shall give without derivation the expression for the eigenvalue Aj(A, {A}) 
of family Ff(X) on the Bethe vector $({A}): 

A/(A,{A})= £ a m (X) N fl c m (A - X k ) . (270) 

m=—s k=l 

The complete list of a m (X) and c m (A) is not important. Suffice to say that 

a m (0) = (271) 
for all rn = —5, — s + 1, ... s — 1 and a s (0) = 1. Further, 

c s (A) = ^— - . (272) 
A + is 

The derivation is based on the relation ( |234| ) which allows to commute the diagonal 
elements of matrix 7/(A) with B(X) from T a (X). From (|270|) we read the momentum 
and energy of quasiparticles 

P (A) = Il„A±| ; (273 ) 
<A) = . (274) 

The relation 

<A) = iA p( A) (275) 

holds; also it is p(X) which enters the LHS of BAE ( |228j ), so we do not need the Lax 
operator L n j(X) to calculate p(X) and e(A). 

The BAE ( |228| ) can be investigated similarly to what was done in case s = 1/2. 
In the limit N — > oo the roots are combined into complexes of type M. The momenta 



of these complexes are distinct of those for spin 1/2. The S- matrices are defined 
by the RHS of BAE and are the same for any spin. One can make the counting of 
roots and get completeness by showing, that full number of states (for which Bethe 
vectors are the highest weights) is equal to (2s + 1)^. 

Finally, we add some comments on the thermodynamic limit in the antiferromag- 
netic case. The ground state $o is given by u s -i/2 = N/2, um — 0, M ^ s — 1/2. 
The excitations correspond to v s -i/2 macroscopic and all other vm finite. They have 
particle interpretation as spin 1/2 particles with the same one particle momentum 
and energy as in the case s = 1/2. However the counting of their states shows, that 
excitations have more degrees of freedom, than just rapidity A and spin e. 

I shall describe the picture of excitations (without derivation) using the language 
of the creation operators. In addition to spin label e and rapidity A this operator is 
supplied by a pair of indeces a, a' assuming integer values from to 2s and subject 
to condition a' — a ± 1. The n-particle excitation is given by a "string" 

Z£'°HAi)3£ ,a3 (A 2 ) . . . Z e a : -^(A n )$ , (276) 

where a = and a n = 0. The counting of Bethe vectors is in exact accord with 
this picture. 

The exchange relation for operators Z employs the S-matrix, which is a tensor 
product of the spin 1/2 S-matrix from section [F| and S-matrix, which acts on the 
indeces a. 

b 



The latter will be denoted by S y x 
tions as follows 



and it enters the exchange rela- 



Z ab (X)Z bc (/2) =J2Z ad (n)Z dc (\)S ( x-n 
d \ 



(277) 



where we suppressed the usual spin variables. They are easily introduced if we write 
the full ^-matrix as 



S = S 1 / 2 ^ 2 (\)®S[ x 



(278) 



The consistency of relations of type ( |277p is based on a star-triangle relation for 

S ( X a d ) 
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The explicit expression for S \ x 



(A) = exp | -i J 



contains factor So (A) given by 



dx sin(Ax) sinh(s + l/2)x 



x cosh(x/2) sinh(s + l)x 



(279) 



In particular 



S x 



a + 1 
a + 1 



a + 2 



S x 



a - 1 



a - 1 



a - 2 



Sb(A) 



(280) 



and 





o + l \ 


(• 


a a J 




tt+1 / 



Sillh (^2( A + + 1 ))) Sil1 2^2 

sin (2^2 ( a + !)) sinh fe( A " z ) 



5 (A) . (28i; 



Other components 

a - 1 
a - 1 



S A 



SI x 



a- 1 
a + 1 



SI x 



a + 1 
a - 1 



(282) 



contain similar trigonometric factor besides Sq(X). 

Thus the antiferromagnetic spin s chain has very remarkable excitations. They 
are particles with rapidity and spin 1/2, but also kinks, relating the "local vacua", 
labeled by a = 0, . . . , 2s. Only transitions between the adjacent vacua are allowed 
in the physical Hilbert space. Some analogy with the Landau-Ginsburg picture in 
the topological field theory is evident, but not explored yet. On this intriguing note 
I finish the general discussion of the XXX S model. 



9 XXX S spin chain. Applications to the physical 

systems 

The appearance of a parameter s in our disposal allows to use it to construct some 
models, going beyond the usual spin chains. Of particular interest is the limit of 
infinite spin s — > oo, combined with the formal continuous limit A — > 0, which can be 
realized in many variants. I shall show, that such representative models of quantum 
field theory as Nonlinear Schroedinger equation (NLS) and S 2 nonlinear a-model can 
be obtained from the XXX S chain in this limit. The first model is nonrelativistic and 
of prime interest in the condensed matter physics, but the second one is interesting 
model of relativistic field theory. 

I begin with the NLS model, and employ the realization of spin variables via the 
complex canonical variables from section [I} 

S n = X*n( 2s ~ XnXn) ] S~ =Xn] S* = XnXn ~ 8 (283) 

and suppose, that s is some complex number. 

The canonical commutation relations for x*, \ n assume the usual form 

[XmiX* n }=5 mn . (284) 

If A is a lattice spacing, Xn and x* n have order A 1//2 in all expressions, where they 
enter in the normal order, i.e. all x* n t° the left of all Xm with the same n. 
The invariant C„ iri+ i of two adjacent spins S n and S n+ i is given by 

2C nin+ i = 2 ^2 S"S" +1 = 2S^S^ +1 + S n S^ + i + S^S n+1 = 

a 

= 2 S 2 - 2s(Xn - Xn+l) (Xn - Xn+l) ~ (X* ~ Xl+lY Xn+lXn- (285) 



The second term in the RHS looks satisfactory, in the formal continuous limit it 
leads to the quadratic form of derivatives of field x{ x )-> X*{ x ) 

X n = A 1/2 x(x), x*n = & 1/2 X*(x), x = nA (286) 

due to the prescription 

Xn+i = A^xOr) + A X \x) + 0(A 2 )) . (287) 

However the last term looks bad and does not lead to the desired expressions 
{x* i x )) 2 (xi x )) 2 characteristic of the NLS model. The remedy is to use equivalent 
variables ip*, ip n on the lattice 

^n=(-lTXn, r n =(-l) n X* n (288) 

and consider ip^,ip n as producing the field ip(x),ip*(x) in the continuous limit as in 
( |286| ) . In the new variables we have 

2C n , n+1 = 2s 2 - 2s(r n + C+iXVv. + VwO + (r n + C+i) 2 </Wn- (289) 

Now the second term in the RHS is bad, but we add and substruct 4s (VO/'n + 
Vn+i^Wi) t° change it into 

2C n , n+ i = 2s 2 - 4s(^> n + ^* n+l ip n+l ) + 

+ 2 S « +1 -V;)(^ + l-^n) + (^ + ^ + l) 2 ^n+l^n , (290) 

which in the continuous limit will contain only good densities n(x) = ho(x) = 
and hi(x) = (if)*) 2 (ijj) 2 . We introduce the operator J via 

J(J+1) = 2s(s + 1) + 2C n ,„ +1 = 

= 4s 2 + 2s - 8sAn(x) + 2sA 3 h (x) + AA 2 h 1 (x) (291) 

and consider the limit s — > oo, A — ► 0, sA = g, where g is a new fixed parameter. 
We see, that the relation ( |291| ) allows to obtain the asymptotics of J in this limit 

a b , 
J = 2s + - + - + ... . 292 

s s z 

Substituting this into the expression 



H n , n+1 (J) = 2^-\nT(z] 



(293) 

Z = l + J 



from section and using Stirling formula for the asymptotics of T(z) and formal 
rule 

A £ = dx (294) 

n J 

we get for hamiltonian the expression 

/fi attice = const --N + ^# NLS + . . . , (295) 

s S J 



where 

N — I i)*i)dx (296) 

# NLS = J [\^'(x)\ 2 + g^f^f] dx . (297) 
This agrees nicely with the dispersion law for the energy of quasiparticles 

m = T^ = -~^ + - , (298) 



A 2 + s 2 s s 3 

if we assume, that the momentum p of NLS particle is connected with the rapidity 
of XXX particle by a simple scale 

p = X/g . (299) 

The state Q plays the role of no particle state of the NLS model. The bound 
states survive in the limit s — > oo when g < 0. For g > the physical hamiltonian 
is usually taken in the form H — p,N, where n is a chemical potential. Here the 
problem of ground state reappears and Dirac sea of particles, created by ip*(x) from 
Q is to be used. The analogous problem for the original magnetic chain consists in 
the inclusion of magnetic field into hamiltonian 

H^H-nS 3 . (300) 

This changes the picture of Dirac sea: only quasiparticles with rapidities, confined 
to some finite interval —B < A < B form the Dirac sea. The excitations are holes 
and quasiparticles with |A| > B. I shall not treat this case in any detail and finish 
the discussion of the NLS model. 

I turn to the second example — the nonlinear a-model with the field n(x) taking 
values in the 2-sphere S 2 . 

The classical field n(x) can be described as a vector n = (ni,n 2 ,n 3 ) subject to 
constraint 

(n,n) = n\ + n\ + n 2 3 = 1 (301) 
for all x. The canonical conjugate variable may be taken in the form 

l = -d nxn, (Z,n) = 0. (302) 

7 

The canonical Poisson brackets 

{l a (x),l?(y)} = e a ^P(x)S(x-y) ; (303) 
{l a (x),nP(y)} = e a ?W(x)5(x-y) ; (304) 
{n a {x)y{y)} = (305) 



and hamiltonian 
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H = jJ^l 2 + , —jdx (306) 
follow from the lagrangian 

C = ^d^nd^n (307) 



and constraint ( |301| ) in a usual way. Here 7 is a coupling constant, which is relevant 
only in quantum case. 

If we regularize the model going to the chain and introducing the variables n k , l k 
as follows, 

l k = Al(x), n k = n(x), x = kA , (308) 
the brackets (|303|) -( j305| ) will assume the ultralocal form 

{Im-i^n} = e 7 ^m^mn ] (309) 

= e a ^nl5 mn ; (310) 
{<,<} = . (311) 

For each lattice site m the phase space is an orbit of the group E(3) of motions of 
R 3 corresponding to the choice of Casimirs 

n 2 = 1; (n,l) = , (312) 

which is also cotangent bundle of S 2 . Corresponding quantum Hilbert space can be 
realized as L 2 (S 2 ). 

To make contact with the spin chains I mention, that this Hilbert space can be 
realized as an infinite spin limit of a Hilbert space of a pair of spin variables of spin s. 
Indeed, comparing the Klebsch-Gordan decomposition, already mentioned above, 

2s 

D S ®D S = J2D 3 (313) 

3=0 

and decomposition of L 2 (S 2 ) 

oo 

L 2 (S 2 ) = J2D j , (314) 

3=0 

we see that 

L 2 (S 2 ) = lim D S ®D S . (315) 

s— >oo 

Thus the a-model variables 1% must be realized through the pair of spin variables 
S2fc-i, SW The most naive way 

Ik = &2k-l + £>2k ; (316) 

n k = 2^(S2k ~ S 2 k-i) (317) 

or inversely 

S 2 k-i = ~ sn k ; (318) 

S 2k = ~l k + sn k (319) 

works. Indeed, from the spin commutation relations we reproduce the first two 
relations (|309|) , (|310|) (in their quantum form) exactly and have 

K, <} = -^e^P8 min (320) 



with RHS vanishing for s — > oo. Furthermore we have 

(l k ,n k ) = 0; As 2 n 2 k + l 2 k = As(s + l) , (321) 

which reproduce ( |312j ) in the limit s — > oo. 

Now the idea is evident: we are to consider the XXX S chain on the lattice of even 
length and take two adjacent points (2k — 1,2k) as one lattice point for a-model. 



Let us see, how the hamiltonian ( |306| ) appears in the formal continuous limit. For 
that we are to estimate the invariants C 2k _i t2k and C 2 k,2k+i- We have 

C 2k -i,2k = -fi - s 2 n 2 k = h 2 k - s(s + 1) (322) 

and 

IS 2 

C 2 k,2k+i = -^hh+i + -(n>klk+i — n k+ih) ~ s rikUk+i (323) 
Now we estimate the invariant in A expansion using the convention ( |308| ) and 

n k+1 = n(x) + An'(x) + ^AV(x) + . . . (324) 

to get 

2C 2fe _ li2fc + 2s(s + 1) = A 2 l(x) 2 ; (325) 
2C 2k , 2k+1 + 2s(s + 1) = A 2 (l{x) 2 -2s(n'(x),l(x)) - s 2 (n" (x) , n(x))) . 

Corresponding operators J 2k ~\ )2k and J 2k , 2k +i have order A 2 so that 

H 2 k-i,2k + H 2k:2k+ i = 2^(1 + J 2k -i, 2k ) + 2^(1 + J 2k<2k+ i)) (326) 
can be easily calculated. Taking the sum over k with our usual understanding 

A ]T = fdx (327) 



and integrating by parts we get for the hamiltonian H a 

H ' - Jwf™ =;/(('- r'i + t'"'' 2 ) * • < 328 > 

which turns into ( |306| ) if we comment, that the map 

/ -> I + an' , n -> n (329) 

for any a is a canonical transformator for the brackets (|303|) -( [3~05"D . The coupling 
constant 7 is connected with spin s as follows 

7 = - . (330) 
s 

The shift ( |329| ) is produced by a topological term, added to the lagrangian ( |307| ), 

•d 

C# = — (cLn x d v n,n)Euv ■ (331) 

07T 



Indeed, such addition changes only the definition of the canonical momenta 



I -»• I + —n' (332) 
Arc v ' 

and thus we can interpret the model obtained from XXX chain as nonlinear cr-model 
with t?-term with 

■8 = 2ns . (333) 

The ^-term with •d = 2im is trivial; in other words •d is defined mod 2tt. Thus we see 
an important difference of integer and half-integer spin s used in our construction. 
The topological term is present only if spin is half-integer. This phenomenon and 
its consequences are discussed in detail by Professor I. Affleck in his Les-Houches 
lectures of 1993 and I can only refer you to the corresponding proceedings. 

Unfortunately the described connection was not yet realized in any real compu- 
tation for the nonlinear a-model. In particular we expect, that the magnetic chain in 
the relevant thermodynamic limit must have an excitation of spin 1 with relativistic 
dispersion law 

p(A) = c sinh A ; (334) 
e(A) = c cosh A (335) 

with parameter c being exponentially small for s — ► oo 

c = e~ s/2 . (336) 

This will lead to the realizations of dimensional transmutation program giving mass 
m via the lattice spacing A and the coupling constant 7 = 2/s 

m = ^-e~ 1/7 . (337) 

Apparently to achieve this we are to understand which portion of the infinite tensor 
product of the spin chain is compatible with the formal continuous limit we just 
described. 

One comment could be relevant to this program. It is natural to combine the 
adjacent Lax operator into product L2k,a(X)L,2h-i, a (X) and perform the change of 
variables (|318|) , ( |319[) there. It turns out, that this new Lax operator has a new 



local vacuum in /^(S ) and the corresponding BAE look like 

/ X k -is X k + i(s + l) \ N/2 = A fc -Xj-i 
\\ k + is X k -i(s + l)J j-^Xk-Xj+i' 

The alternating value of spin necessary to maintain our conventions ( |3Q8| ) is manifest 
here. However the proper choice of solutions to these equations is not done yet. I 
leave it as a challenge and stop here the discussion of nonlinear cx-model. 

I also finish considerations of XXX chains (with one revisit in section |12]) and 
turn to the XXZ chains. 



10 XXZ model 

As was told above, this model is a deformation of XXX model with one new param- 
eter. We shall denote this parameter by q or 7 with q = e n . The deformation uses 
the q- analogues of usual number 



q x — q x sin 7X tt / x + n7r 7 1 



1 r, n-X 



q- q 1 sin 7 ^ V 1 + rm^y 1 



which effectively change the complex plane to a strip via the multiplicative averaging 
in (|339|) . We begin with the construction of the Lax operator of XXZ model using the 
matrix analogue of this averaging. For classical spin variables the formal averaging 

00 

^f(A)= II ^(A + ifory- 1 ) (340) 

k=— 00 

can be evaluated to lead to the expression 

« Z (A) = -^f Sinh7(A ^ } . w7%3,) > ^ 341 ) 
sin 7 \ SU17S+ sinh7(A - zS£) J 

where are some function of the original spin variables S„, of XXX model. 

Taking this as heuristic consideration we shall look for the Lax operator L na (X) in 
the form 

(xq^ n — x~^q~'~* n (g — q ~~^)*S'~ \ 
1 -i\ c+ -s 3 -i™5 3 ' (342) 
(q — q )b£ xq J ™ — x q^ n J 

using multiplicative spectral parameter 

x = q~ iX (343) 
and quantum operator g s ™, S^. We shall check, that L n>a {x) satisfy the FCR 

R ai ,a 2 ( x /y) L n ,ai(x)Ln,a 2 (y) ^ L nta2 (y)L nm (x)R aija2 (x/y) , (344) 
where R is a g-deformed analogue of the R- matrix from section |] (see (f73|) ) 



R 



(345) 



fa \ 
b c 
c b 
V a J 

with 

a = qx — q~ x" 1 ; b = x — x~ 1 ; c = q — q~ 1 . (346) 
To check FCR (|344|) it is convenient to twist it a little, introducing 

L n , a (x) = Q(x)L n , a (x)Q-\x) , (347) 
R ai , a2 (x) = Q{x)®Q{y) R ai , a2 {x/y)Q-\x)®Q-\y) , (348) 



where Q(x) is a matrix in the auxiliary space 

,1/2 

x- 1 ' 2 



Q(x)= X T Jm )■ (349) 



It is clear, that FCR is true for original L n a (x) if it holds for L n ^ a (x). 
Now observe, that L n>a (x) and R ai) a 2 { x ) have a simple x-dependence 

L = xL + - x- x L_ , R = xR+ - x~ l R. , (350) 

where we dropped indeces n, a±, a 2 ,~, and matrices L + , L_, R + , i?_ are given by 







(q - q-^S- 
q-S 3 



Ra 



1 q — q 1 
1 



q- s * 
-{q-q- 1 )S + q s " 

\ 



(351) 



(352) 



q J 



( q- 1 



1 

iq-q- 1 ) 1 



q- 1 J 



It is clear, that 
and 



R- = Pi?; 1 ? 



(353) 

(354) 
(355) 



R + -R- = (q- q- 1 ) P , 

where P is a permutation matrix. 

Separation of spectral parameters in FCR shows, that it is true, if the following 
7 relations hold: 

RL\L 2 ± = L\L\R , (356) 
where R can be R + and R_ and labels 1 and 2 substitute a\ and a 2 ; furthermore 



1 T2 



R+L^L 

r_i}l\ 



and 



R+L^ L 2 , 



r_l\l 2 



L_L + R^ 
L 2 L l R_ 



LiLR-i 



L 2 L\R_. 



(357) 
(358) 

(359) 

Only three of them are independent and we take as such two of the relation ( |356| ) 
for L + and L_ separately and relation ( 357 ). The two other relations in ( |356|) and 
relation ([35 8|) are easily reduced to the chosen ones if one applies the permutation 
taking into account, that 

L 1 = PL 2 P . (360) 

Finally relation (|359|) is checked if one uses the property (|355|) to substitute unwanted 
R+ by R- and vice versa. 

Now it is easy to check, that the basic relations, say 



R+L^L 2 , 
R+L^L 2 ^ 

r+l\l 2 



lP^l.r^ 
l 2 _l\r^ 

L 2 _L\R a 



(361) 
(362) 
(363) 



are satisfied, if q s3 , S + and S satisfy the commutation relations of the g-deformed 
si (2) from section [| This finishes the proof of the FCR ( |344| ) . It is worth to mention, 



that the XXX S Lax operator (|218|) is obtained from (|342 ) in the limit q — > 1, so 
FCR is proved also for it. 

Another comment is, that we can consider the relations (|361|) - (|363|) together 
with the structure (|351|) as alternative definition of the g-deformed sl q {2) algebra, 
which will prove to be convenient in what follows. 

If we renormalize R± by 

R+ q~ 1/2 R+, R- -> q' 1/2 R- (364) 

we see that they take form ( |351| ) where q s3 and realize the 2-dimensional repre- 
sentation 

deforming the Pauli matrices. 

The i?-matrix and Lax operators above can be considered as representation of the 
universal 7£-matrix, which in this case corresponds to the g-deformed afline algebra 
sl q (2). 

In some sense this algebra (see a lot about it in Professor Miwa lectures) is more 
natural object than Yangian, which is a special contraction of it when q — > 1 (or 7 
tends to 0) in conjunction with a renormalization of the grading parameter A. 

The g-deformed si (2) has finite dimensional representations, analogous to those 
in the nondeformed case; their dimension is 2s + 1, where spin s is integer or half- 
integer. Besides it has other interesting representations, called cyclic and already 
briefly mentioned in section § these representations have no limit, when q — > 1. 

The deformed representations in C 2s+1 are of highest weight type, so that there 
is a state uj such, that 

S + co = 0, q s3 uo = q s uo . (366) 

This allows to repeat all constructions of Algebraic Bethe Anzatz and the BAE 
assume the form 

sinhfAfc + isj) _ smh(\ k - Xj + ij) / 367 \ 
sinh(A fc - is'y) sinh(A fc - \ - ij) 



where we returned to the additive spectral parameter A to make ( 367|) look rather 
similar to BAE (ggg). 



The investigation of these equations is to large extent similar to what was done 
before in XXX case. One encounters complexes of roots, one can estimate the 
number of Bethe vectors etc. A new feature is a role, played by the arithmetic 
nature of 7. For instance if 7/7? is rational then q is a root of unity and all specifics 
of sl q (2) for such a case is to be taken into account. We have no time to discuss it 
in any detail. 

I turn now to the problem of constructing the local hamiltonian. For that we 
need to control better the invariant of the pair of spins. 
It is convenient to use the matrix 



L = L + LZ 



(368) 



to describe the deformed spins. The matrix L_ Y is given by 

£: '=U-"V A)- (369) 

The contractions q — > 1 is especially transparent here with expansion 

L = I + 2i 1 (^ S s l ) +... (370) 

for 7 — > 0, where S" 3 , 5 are nondeformed spin variables. 

The entries of matrix L satisfy the relations which can be cast in the matrix 
form 

L 1 RZ 1 L 2 R_ = R^-tfR+L 1 . (371) 
There exists the generalized trace (so called g-trace) 

ti q A = ti(DA) (372) 

with 

such that for any matrix A we have 

tr 2 q R- 1 A 2 R = ti q l 1 A , (374) 

where we calculate tr 9 only over the second factor in the tensor product of two 
auxiliary spaces. 

From these properties we see immediately that 

[tr g L, L\ = , (375) 

i.e., tr g L plays the role of Casimir and we shall denote it as C. 
The explicit calculation gives 

C=(q + q- 1 ) (q 2s3 + q~ 2s3 ) + (q - g" 1 )' (S + S~ + SS + ) , (376) 

or 

C = q 2s3+1 + q- 2S ' i - 1 + 2(q-q- 1 ) 2 S-S + , (377) 

so that in the irreducible representations of spin s the eigenvalues of C are given by 
q2s+i _|_ q-2s-i fj^jg prompts the introduction of the operator J such that 

C = q 2J+1 + q- 2J -\ (378) 

which is an analogue of J from the usual sl{2). 

Armed with this knowledge we can attack the problem of local hamiltonian in a 
way analogous to that of section || 

We shall take two spin operators and describe them by the corresponding L- 
operators, which we shall denote by L± and M±. Let L(x) and M(y) be the corre- 
sponding Lax operator, taken in the form fl350|) . The FCR we shall use to find the 
representation for the R s ' s (x) Lax operator R(x) looks as follows 

L (1/x) M (l/j/) R (x/y) = R (x/y) M (1/y) L (1/x) . (379) 



Putting R (x) into the form 

R{x) = Pr(x) (380) 

we have instead 

M (1/x) L (1/y) r (x/y) = r (x/y) M (1/y) L (1/x) . (381) 

Altogether ( |381| ) comprises four matrix equations. Taking r to be a function of 
invariant of spin S and T we trivially satisfy two of them. One of remaining looks 
as follows 

(x^M+L- + a;M_L+) r(x) = r(x) (xM+L„ + x^M-L+j . (382) 
We shall take lower left element of this equation 

(x- 1 q- T3 S + + xq s ' A T + ) r(x) = r(x) (xg" T3 S + + x~ l q s *T + ) (383) 

and consider it in the subspace 

g T3 S + +TV s3 = , (384) 

where the invariant takes the form ( [378 ) with 

J = T 3 + S 3 _ ( 3g5 ) 

We shall look for r(x) as a function of q 2J+l . Our equation looks as follows 
(x- 1 q- J - xq J ) q s3 S + r(x, q 2J+1 ) = r(x, q 2J+l ) (xq- J - x^q^ q s " S + . (386) 
Using the commutation relation 

S+q 2J+1 = q- 2 q 2J+1 S + (387) 
we transform the equation ([386 ) into the functional equation 

r(qw) 1 — x 2 w 



r(q 1 w) x 2 — w 

where 



(38^ 



w = q 2J+l . (389) 

We shall not discuss this equation here; suffice to say, that it is a g-deformation of 
equation ( |253| ) and its solution is given by some g-deformation of T-function (known 
also as g-exponent and exp of quantum dilogarithm) . This finishes the general 
discussion of the XXZ model. 



Figure 1: Discrete space time and initial saw 



11 Inhomogeneous chains and discrete time shift 

Here I shall present some development of the general scheme of Algebraic Bethe 
Anzatz which will allow to include more dynamical models under its spell. Until 
now two main observables — momentum and energy were treated differently. We 
used discretized space, but continuous time; thus we introduced a finite space shift 
U and infinitesimal generator of time shift H . To make our consideration more 
manifest invariant it is natural to discretize also time. Exactly this will be done here. 
We used lattice as a way of regularization and always had in mind corresponding 
continuous limit. But some people take discrete space-time more seriously as an 
inevitable consequence of gravity. I shall not open this discussion here sticking to 
the formal aspects only. 

The shift operator U was obtained above as a trace of the monodromy at some 
distinguished value of spectral parameter. We need now two such distinguished 
values. The way to achieve this goal is to consider a specially inhomogenous chain 
with the spectral parameter taking alternating values A ± uj for some fixed uj. This 
simple idea indeed works as we shall see momentarily. 

Let L n j(X) be a fundamental Lax operator with the same quantum space h n 
and auxiliary space Vf. The monodromy of the inhomogeneous chain is given by 

T f (X, uj) = L 2NJ (X + w)L 2 N-i,f(* -<*>)... L 2,/(A + uj)L xj {\ - uj). (390) 

We want to argue, that 

U+ = tT f T f (uj,uj), U- = tT f T f (-u,uj) (391) 

play the role of shifts in the characteristic (light-like) directions in the discrete space- 
time. This space-time is natural to draw as shown in Figure [3], where the monodromy 
is taken as a product along the initial saw, depicted by a fat line. 

The space runs from east to west and time from south to north. The L2 n ,/(A + c<j) 
is transport along NW direction and L2 n +i,/(A — uj) is the same in SW direction. 

To see the shift properties of operators U± we consider even more inhomogeneous 
monodromy inserting into the string of Lax operators L n j{\ ± uj) the Lax operator 
Lf^(fj, — A), treating the space a as quantum and / as auxiliary. This admixture is 
put between L 2n j(X + uj) and L 2n -i,/(A — uj). The new monodromy looks as follows 

T/(A, u\a, n, fx) = L 2 7V,/(A + uj)L 2 n-ij(^ - uj) . . . (392) 
• • • L 2n ,f{^ + uj)Lj^ a (n - A)L 2 „-i,/(A -uj)... L 2 ,f{\ + uj)Lxj(X - uj). 




Figure 2: Zero curvature equation 

Due to the similar FCR for all entries here this monodromy also satisfies FCR and 
its trace 

F(X,u\a, n, /x) = trjTy(A, u\a, n, /j) (393) 

gives a commuting family as function of A. 

Now let us look at F(X, u\a, n, fi) at distinguished points A = ±u. We use the 
normalization 

L nJ (0)=P nJ , (394) 
where P is a permutation of h n and Vf, which coincide as spaces. We have 

F(u,u\a, n, A) = tr/ (L 2NJ (2uj)P 2N ^ 1 j . . . 

. . . L 2nJ (2cu)Lj'(\ - u)P 2n ^j . . . Uj^P.j) . (395) 
Note, that we changed parameter \x to A. Due to relation 

Lf a {\ - u)P 2n „u = P 2 n-i,/£ 2 t 1>a (A - u) (396) 
we can bring LJ^(X — to) to extreme right to get 

F {u,u\a, n, A) = U + L^_ 1>a {X - u). (397) 

Analogously 

F {-u,u\a, n, A) = L^ a (\ + u)U-. (398) 
The commutativity of F in the LHS of ( |397| ) and ( |398| ) leads to the equation 

U + L 2 l ha (\ - u)L^{X + = L^ a (X + u)U-U + L^_ ha (\ - w ) , (399) 

which we can rewrite due to commutativity of U + and U- in the form 

£ 2 n,a(A + u)L 2n _ lja (\ -u) = *7_L 2n _ M (A - t^CT 1 £/+%n,a(A + w)U+. (400) 

This equation has natural interpretation as a zero curvature condition for the trans- 
port around the elementary plaquette in our space-time. On Figure we use nota- 
tions 

^2n-i,/(A - u) = U^L 2n ^j{\-uj)UZ l ; (401) 
L 2nJ (\ + u) = U- l L 2nJ (\ + cu)U + . (402) 



Thus we see, that U+ is a shift in NE direction and U_ 1 is shift in NW direction. 
In terms of shifts in N (time) and W (space) directions e~ lH and e~ tP we have 

U+ = e~^ H ^ 2 , U. = e WV 2 . (403) 

This is our main assertion. Now I shall describe a more explicit expressions for U±. 
Using our usual Ansatz 

L n>f {\) = VnJnjW (404) 

we get 

U+ = tT f (P2N,fl2N,f( 2 ^) P 2N-lJ • • • V 2J l 2J {2uj)Y 1 j) (405) 

and bringing all l n j(2u) to right we get 

U+ = V'[[l 2nM - 1 (2u) , (406) 

n 

where V~ x is a shift n — > n + 1. 
Analogously we have 

U- =II / 2n,2n-l(-2^)\/ , (407) 

n 

kn,2n-l{-\) =^Vl(A) • (408) 

From identification ( |403| ) we get 

e iP = U + U_ = V 2 ; (409) 

= U+UI 1 = VY[l 2n ,2n-l(2uj)V- 1 lll2n,2n-l(2uj) = 

= n / 2n+l,2n(2^)n^2n,2n-l(2cu) . (410) 

We see, that the physical space shift is shift n — > n + 2; in other words two lattice 
sites of our space lattice constitute one physical site. We already have seen such 
trick in the discussion of nonlinear a-model. 

The derivation of the BAE for the inhomogeneous chain does not differ from that 
given above. The equations look like 



' a(Xj + u)a(Xj — uS) \ N 
v 8(Xj + cu)5(A j - u) J 



n^-A,) , (411) 

Mi 



where a(X) and 5(A) are local eigenvalues and factor S(X) in the RHS is a quasipar- 
ticle phase factor. From this we read the quasiparticle momentum and energy: 

ip a(X + u)a(X - uj) ih 5(X + ij)a(X - u) 

5(X + Lu)5(X-Lu) , a(X + u)5(X-uj) ' { ' 



These expressions nicely turn into our previous formulas in the limit uo — - > 0. This 
finishes our description of a general scheme and I turn to examples. 



12 Examples of dynamical models in discrete 
space-time 



I shall present two examples for the general scheme of the previous section. The 
first is associated with XXX S spin chain. The BAE 

V ( A i + u ~ is )( x i -u-is)J k}j - A fc - i 

can be investigated as above. The spectrum around the antiferromagnetic state, 
already described in section |8], leads to the dispersion law for the excitations 

p(\) = arctg sinh 7r(A + uj) + arctg sinh 7r(A — uj) ; (414) 

e ( A ) = z n t + u n (415) 

cosh 7r(A + uj) cosh7r(A — uj) 
In the limit uj — > oo we get the relativistic one particle spectrum 

p(X) = m sinli7rA , e(A) = m cosli7rA , (416) 
where m is obtained via the "dimensional transmutation" 



1 

A" 



(417) 



if we introduce lattice spacing A to anticipate the continuous limit. 

Let us concentrate on the spin of the excitations. The results of section || are 
rather complicated. However they drastically simplify in the limit s — ► oo. Indeed 
the restriction a <2s for kink labels disappears and the sequence 0, a 1; . . . , a 2n ~i, 
can be considered as parametrizing a particular singlet, entering the representation 
of sl(2) in n 2n ®C 2 . For n — 1 (two particle state) we have just one such state 
(0, 1, 0), corresponding to a singlet in C 2 ® C 2 . The phase factor, corresponding to 
the triplet state in spin space and singlet in the kink space is given by the s —>■ oo 
limit of (|28 1|) and (|279|) and is rather simple: 



S(X) = S t (\)^:S t (\) , (418) 

/A | % 



where «St(A) is given by ( |211|) . The last expression is nothing, but the phase factor 
for a (triplet) ® (singlet) state with respect to two independent s/(2) groups with 
S'-matrix, given by a tensor product of two 5*- matrices of type S 1 ^ 2 ' 1 ^ 2 . 
Now we mention, that 

sl(2) ® sl(2) = o(4) (419) 

with C 2 <8> C 2 being a vector representation for it. In other words, we can interpret 
the excitations in our inhomogeneous XXX S model in the limit s —>■ oo as vector 
particles corresponding to the vector representation of o(4). 

This particle content coincides exactly with what is believed to be true in the 
S 3 nonlinear a-model (or the s/(2) principal chiral model). The inhomogeneous 
XXX S spin model just realizes one particular sector of sl(2) chiral model in the limit 



s — > oo; only particles in singlet state with respect to the right spin appear. However 
this is enough to characterize the full S-matrix. 

There are more reasons to justify this correspondence, in particular the exact 
calculation of the /^-function of the renormalization group is possible via the BAE. 
I do not have time to discuss it here. 

The second example is the Sine-Gordon model. This dynamical system was 
instrumental for the developing the ABA by Sklyanin, Takhtajan and me in the end 
of 70-ties. From that time there were developed quite a few approaches to investigate 
it via BAE. Here I present an approach based on alternating inhomogeneous chain, 
developed by A. Volkov and me at 1992. 

As the dynamical variables on an alternating lattice I shall use the Weyl variables 
u n , v n with the exchange relations 

u n v n = qv n u n , q = e n . (420) 

The auxiliary Lax operator 



belongs to the XXZ family with FCR (|344j) . I cannot help mentioning, that in 



simplicity formula ( f421| ) beats even the expression for the Lax operator of XXX 
model fl32]). 

The local product of two Lax operators L2 n ,a(xK,)L2 n -i,a{xK~ 1 ) corresponds to 
transport along a physical lattice site. We have explicitly the matrix 



—x 



U 2n U2n-l ~ X 2 V 2n V 2 n-l X ( K 1 U 2n V 2n -l + KV 2nU 2 n-l) 



(422) 



We can identify this matrix with the XXZ Lax operator of type ( j342|) with the spin 
variables q si , given by (|22|), ( p3|) if we impose the constraint 

W n = U 2n U 2n -iV 2n V2n-l = 1 (423) 

at each physical lattice site. This constraint reduces the number of local degrees of 
freedom from two to one, as it must be. 
Indeed, let us put 

e inn = u 2n v 2n - x ; (424) 
e iVn = u 2n u 2n -x ■ (425) 



Then ( |422[) turns to ( |342| ) after division by x and similarity transform with matrix 

D = ( ^ K % ) , (426) 



if we put 



2 2 
771 = K 



(427) 



After this identification it is natural to assume that the BAE look like (|411 ). The 
derivation is nontrivial, because the representation for the g-deformed spin variables 
is not of the highest weight. However it can be done and the local factors a(X) and 
5(X) correspond to spin —1/2. Returning to the additive variables x = e x , k = 
we have the BAE in the form 

/ sinh (A,- + u + f ) sinh (A, - u + g ) \ N/ * = i sinh(A, - A, - ij) 
\sinh (\j + u) - f ) sinh (Aj - u - f ) / " fc= i sinh(A j - A fc + 17) ' 

Now I turn to the description of the time shift operator. For that we are to find 
the fundamental Lax L n j(x) operator, corresponding to the auxiliary L n>a (x) from 
([42 1|) . We write the equation 

L 2 (1/x) U (l/y) r (x/y) = r (x/y) U (l/y) L 2 (1/x) (429) 

with the natural brief notations 

Li(x) = L nua (x), L 2 (x) = L n ^ a (x) . (430) 

The solution will be looked for as a function of the variable 

w = u 2 u\v 2 x v\ . (431) 



The diagonal matrix elements of equations ( |429| ) are trivially satisfied. The right 
upper matrix element of the equation ([42 9|) looks as follows 

(xu 2 v i + v 2 Ui 1 )r(x, w) = r(x, w)(u 2 vi + xv-iu^ 1 ) (432) 

and can be reduced to the functional equation 

r(x, qw) xw + 1 



r(x, q x w) x + w 



(433) 



by using the exchange relations ( |420|) . We have already seen this equation in section 
[T0| . Indeed, for the representation (p2|), ( p3| ) we can check, that w coincides with 

g2J+l 

Denoting 

w 2n = u 2n u 2n „iv^v 2n „i ; (434) 
w 2n+ i = u 2n+ iu 2n v 2 n +1 v 2n , (435) 

we have the expression for the time evolution operator 



e 



- Y\r(K\w even )\{r(K\w oAA ) . (436) 



This expression is simple enough to allow us to write down the equation of motion 
explicitly. It is convenient to introduce a new set of variables ip w such that w n are 
their "multiplicative derivatives" 

w n = . (437) 

Vn-l 




Figure 3: Elementary plaquette 

Variables i]) m and ip n commute for n — m even, so there is no problem in order of 
factors in ( ]437| ) . The exchange relations among ip n are nonlocal, but it is not of 
concern to us; it suffices to say, that given ip n does not commute with only one w, 
namely w n , and for this pair there is the Weyl relation 

ip n w n = q 2 w n ip n . (438) 

We can take the set {Vw n , w evea } or {^odd, w dd} as a set of independent Weyl 
variables in the physical Hilbert space. I skip here the subtle question of boundary 
conditions for integrating (|437|) , which is settled in the original literature. 

Equations of motion are especially simple in if) variables. Let ip n be a variable 
ip n once displaced in time 

^ n = e lH iP n e~ lH . (439) 

Take for defmiteness ip Q dd, i-e. ip 2n +i- Only one factor, containing w 2n +i in ( |436| ), 
does not commute with it, so we have using the functional equation (|433|) 



^2n+l = r 1 (w 2n +l) 4>2n+l r(lU 2n +l) = V>2n+l r 1 (q Vin+lM^n+l) = 

= ^2n+l 2 ^ Zi = ^n+l—T, ~ 2 • (440) 

This is our equation of motion. For ip even the derivation is exactly the same. On the 
discrete space-time drawn on Figure the equation (|440|) connects the ip attached 



to an elementary plaquette shown in Figure |^ and can be rewritten as 

^ N = ^s^Ti — 7 2 / ( 441 ) 

or 

(q^ipN^Pw ~ ^s^Pe) = K 2 (q~ 1 i's^w ~ ^n^e) ■ (442) 
The last simple quadratic equation looks quite appealing. 

Let us show, how the Sine-Gordon equation appears in the formal continuous 
limit. We also confine ourselves to classical commuting variables. 

The naive prescription that ip n define a continuous function of space-time vari- 
ables does not work. We have already seen a similar situation in the case of NLS 
model. The way out is to use variables 

X = { (443) 



alternatively on each second SE characteristic line of our lattice. Then equations 
of motion (|441|) transform into 

- l K 2 q~ l XwXE + 1 
Xn = Xs -=T t~2 444 

q 1 xwXe + k 2 

and for large k 2 and classical limit q = 1 we have 

^ = l + -J— XEXw]+... • (445) 

XeXw « 2 VXbXw / 

We put here 

X = e iip (446) 

and consider as defining a smooth function of x,t in the continuous limit. Then 
evidently 

XNXS = exp { i—(<ptt - (fxx) + ■■■} (447) 
XeXw { 2 J 

and introducing the rescaling 

= m 2 A 2 (448) 



we get from (|445|) the Sine-Gordon equation 

- ^ + 2m 2 sin 2y? = . (449) 

In quantum version the scaling (|448 ) is to be modified, accounting to the mass 
renormalization. The investigation of BAE Q428 ) allows to describe this renormal- 
ization exactly, adding to the nonperturbative intuition in Quantum Field Theory. 
However I cannot speak about it here and stop technical developments to come to 
some conclusions. 



13 Conclusions and perspectives 

The first and the last formulas in this text being the same, our exposition closed 
the circle. We have learned the technique of Algebraic Bethe Ansatz for solving 
integrable models and have shown, how it works in detail on the simplest example 
of spin 1/2 XXX magnetic chain. Several other models were treated more superfi- 
cially, only the specific details were given. Several parameters, appearing in these 
generalizations: spin s, anisotropy parameter 7, shift u in the alternating chain, 
allow to include in our treatment most known examples of soliton theory, including 
relativistic model of Quantum Field Theory. Thus the spin chains showed their 
great universality. If we add, that the fashionable models of conformal field theory 
also can be included as particular limits, then we could claim, that we deal with the 
classifying object in the theory of integrable models. 

We treated here only rank 1 case associated with the Lie algebra sl(2). Gen- 
eralization to the higher rank simple Lie algebras is possible and many results are 
already known. In particular, one can use a Cartan subalgebra label to introduce 
one extra dimension for space. 



We did not make a stress on the subject of quantum groups or quantum sym- 
metries, as this was discussed in other lectures at the School. However I suppose, 
that I made it clear, how quantum groups were created naturally inside of ABA (of 
course, their appearance in the C* algebra approach of Professor Woronowitz was 
completely independent). 

Discrete space-time approach to quantum integrable models seems to me most 
promising and elegant. It has already allowed to write explicitely the quantum 
equation of motion. I cannot help feeling that it is just the beginning of analysis, 
including complex analysis, on discrete manifolds. 

From pure mathematical point of view many of reasonings in my lectures were 
only heuristic. The statements were "claims" rather than "propositions" not to 
mention "theorems" . I am not a supporter of absolute rigor in mathematical physics, 
but would be happy, if some of my claims become better justified. Even more 
important is the fact, that new mathematical constructions and results could appear 
in course of this justification. Especially promising are problems of infinite tensor 
products in antiferromagnetic case and their connections with infinite-dimensional 
Lie algebras and analysis on the discrete space-time. I leave these challenges to my 
listeners. 



14 Comments on the literature on BAE 

Here I will add some comments which will serve as a guide to the literature on Bethe 
Ansatz. Only seminal (from my own point of view) papers will be mentioned. The 
recent monograph |l[ contains quite vast list of references. Another monograph || 
reflects the developments around Bethe Ansatz mostly prior to the advent of ABA. 

BAE were first written by H. Bethe in [[| for the spin 1/2 XXX model. He 
did not use the uniformizing rapidity variables, which apparently appeared first in 
a paper by Takahashi j|] . First integral equation for the distribution of roots was 
derived by L. Hulthen []5J. My exposition of the spin 1/2 XXX model follows closely 
my paper with Takhtajan ||. The statement, that the spin of a spin wave is 1/2 
was announced by us in 0. The scattering theory for the excitations over the 
ferromagnetic vacuum was developed by Babbit and Thomas ||. 

General formalism of the ABA was worked out by Sklyanin, Takhtajan and me 
in U on the example of the Sine-Gordon model. The state of art at the end of 



70-ties was described in my survey [10 



The g-deformed si (2) algebra with defining relations flT2]), flT3| ) was introduced by 
Kulish and Reshetikhin in 1981 in [0]] in connection with the higher spin XXZ model. 
More general investigation of the quadratic algebras was performed by Sklyanin fT2" 



who pointed out the connection with the Hopf algebras |fL3| . The proper algebraic 
understanding of general structures, appearing here, was formulated by Jimbo [0 
and in the most perfect form by Drinfeld fI5f . The role of the fundamental Lax 
operator was underlined by Tarasov, Takhtajan and me in [Q. There we followed 
the ideas of fusion, introduced by Kulish, Reshetikhin and Sklyanin in |nj, where 
the formula (|254|) was derived. 



Investigation of BAE for higher spin XXX model was done by Takhtajan JT8] |. 



see also Babujan |T9|)- The kink interpretation of the excitations was developed 



by Reshetikhin in 
chemical potential was done by Lieb and Liniger f21 



The thorough investigation of BAE for NLS model with 



The interpretation of the NLS and SG models as XXX and XXZ chains 

[0 



cor- 
respondingly, was done by Izergin and Korepin |p2|| . Korepin was also the first 
to propose the formula of type ( |211| ) for the matrix element of .S-matrix for the 
excitations above Dirac sea 



23 



The functional equation for the fundamental Lax operator (|388|) for XXZ model 
has a natural place in the theory of sl q {2) algebra, see e.g. |23[. 

The idea of imbedding the S 2 a-model into XXX chain belongs to Haldane [25 



and Takhtajan and me |[26|| . More detailed development was done by Affleck, who 
in particularly stressed the role of the topological $-term |T7f| . The BAE ( |338|) was 
derived by Bytsko (unpublished). 

The idea of the usefulness alternating inhomogenous chain to treat the relativistic 
models was underlined by Reshetikhin and me |28| in our treatment of the principal 
sl(2) chiral model. It was developed further by Destri and De-Vega p9f . The zero 



31 



curvature interpretation was given by Volkov and me in [30]. 

The quantum equations of motion ( [442| ) for the SG model appeared in 
they coincide in the classical limit (q = 1) with Hirota equations ||32|| . In terms of 
variables f A , connected with ip A by 



Ia - f. 



B 



B 



in the form 



(/AT — fE){fw — fs) 



(450) 



(451) 



U'n ~ fw)(fE — fs) 

as was commented by Volkov ||33|| . This makes the appearance of discrete complex 
analysis indispensable. The same equations were recently advertised by Capel et al. 
[0 . Another line of thought on discrete geometry and discrete SG equation belongs 
to Pinkal et al. [35|. Quadratic algebras appear now in many instances and guises. 
I mention the discretized affine algebra, introduced by Semenov-Tjan-Shansky and 
discussed in some detail in JJfJ. 

I finish by reference to my previous lecture course |37]] where many aspects of 
these lectures were treated, and to lecture notes of my Schladming lectures on new 
applications of BAE to Hofstadter model from condensed matter physics, to high 
energy QCD and to Liouville model of Conformal Fiel Theory. 

In my lectures I did not refer at all to the parallel development in classical 
statistic physics. The general equivalence exists between 1 + 1 dimensional quantum 
dynamical systems and 2 dimensional models of classical statistical physics. Inte- 
grable models of the former domain correspond to the exactly soluble models in the 
latter. A lot of well known results here are connected with the names of Onsager, 
C. N. Yang, Lieb, Baxter and others. I can refer to the monograph of Baxter |38], 
where this theme is displaced in great detail. 

Another connection with statistical physics is the use of the integrable hamil- 
tonians to define the corresponding Hibbs state z~ 1 e~^ H+ ' lN . This was pioneered 
in paper of C. P. Yang and C. N. Yang ||39|| , which led to important development, 
called the Thermodynamic Bethe Ansatz. Already mentioned paper [^] of Babujan 
was an important step in formulating this technical development. 



Returning to quantum field theory interpretation I must say, that until now ABA 
has given in the way to investigate mass spectrum and ^-matrix. More detailed off 
shell characteristics of systems under consideration are very scarce. The main success 
was achieved by Smirnov in the discussion of formf actors of local operators [3D]. 



Another line of thought, originated by Korepin, is partly described in This 
topic is the most important from the point of view of mathematical physics and 
physical applications. 
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